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METHOD FOR CREATING POLYNUCLEOTIDE AND 
POLYPEPTIDE SEQUENCES 

5 

CROSS-REFERENCES TO RELATED APPLICATIONS 
This application derives priority from USSN 60/067908, filed 
December 8. 1997. which is incorporated by reference in its entirety for all purposes. 

0 TECHNICAL FIELD 

The invention resides in the technical field of genetics, and more 
specifically, forced molecular evolution of polynucleotides to acquire desired properties. 

BACKGROUND 

5 A variety of approaches, including rational design and directed evolution, 

have been used to optimize protein functions (1,2). The choice of approach for a given 
optimization problem depends, in part, on the degree of understanding of the relationships 
between sequence, structure and function. Rational redesign typically requires extensive 
knowledge of a structure-function relationship. Directed evolution requires little or no 

-0 specific knowledge about structure-function relationship; rather, the essential features is a 
means to evaluate the function to be optimized. Directed evolution involves the 
generation of libraries of mutant molecules followed by selection or screening for the 
desired function. Gene products which show improvement with respect to the desired 
property or set of properties are identified by selection or screening. The gene(s) 

25 encoding those products can be subjected to further cycles of the process in order to 

accumulate beneficial mutations. This evolution can involve few or many generations, 
depending on how far one wishes to progress and the effects of mutations typically 
observed in each generation. Such approaches have been used to create novel functional 
nucleic acids (3. 4), peptides and other small molecules (3 », antibodies (3). as well as 

30 enzymes and other proteins (5, 6, 7). These procedures are fairly tolerant to inaccuracies 
and noise in the function evaluation ( 7). 
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Several publications have discussed the role of gene recombination in 
directed evolution (see WO 97/07205. WO 98/42727, US 5807723. US 5,721.367, 
US 5,776,744 and WO 98/41645 US 5.81 1,238, WO 98/41622, WO 98/41623. and 
US 5.093.257). 

5 A PCR-based group of recombination methods consists of DNA shuffling 

[5, 6], staggered extension process [89. 90] and random-priming recombination [87]. 
Such methods typically involve synthesis of significant amounts of DNA during 
assembly/recombination step and subsequent amplification of the final products and the 
efficiency of amplification decreases with gene size increase. 
1 0 Yeast cells, which possess an active system for homologous 

recombination, have been used for in vivo recombination. Cells transformed with a 
vector and partially overlapping inserts efficiently join the inserts together in the regions 
of homology and restore a functional, covalently-closed plasmid [91]. This method does 
not require PCR amplification at any stage of recombination and therefore is free from the 
15 size considerations inherent in this method. However, the number of crossovers 

introduced in one recombination event is limited by the efficiency of transformation of 
one cell with multiple inserts. Other in vivo recombination methods entail recombination 
between two parental genes cloned on the same plasmid in a tandem orientation. One 
method relies on homologous recombination machinery of bacterial cells to produce 
20 chimeric genes [92]. A first gene in the tandem provides the N -terminal part of the target 
protein, and a second provides the C-terminal part. However, only one crossover can be 
generated by this approach. Another in vivo recombination method uses the same tandem 
organization of substrates in a vector [93]. Before transformation into E. coli cells, 
plasmids are linearized by endonuclease digestion between the parental sequences. 
25 Recombination is performed in vivo by the enzymes responsible for double-strand break 
repair. The ends of linear molecules are degraded by a 5 , ->3' exonuclease activity, 
followed by annealing of complementary single-strand 3' ends and restoration of the 
double-strand plasmid [94]. This method has similar advantages and disadvantages of 
tandem recombination on circular plasmid. 

30 

SUMMARY OF THE INVENTION 
The invention provides methods for evolving a polynucleotide toward 
acquisition of a desired property. Such methods entail incubating a population of 
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parental polynucleotide variants under conditions to generate annealed polynucleotides 
comprises heterodupiexes. The heteroduplexes are then exposed to a cellular DNA 
repair system to convert the heterodupiexes to parental polynucleotide variants or 
recombined polynucleotide variants. The resulting polynucleotides are then screened or 
5 selected for the desired property. 

In some methods, the heterodupiexes are exposed to a DNA repair system 
in vitro. A suitable repair system can be prepared in the form of cellular extracts. 

In other methods, the products of annealing including heterodupiexes are 
introduced into host cells. The heterodupiexes are thus exposed to the host cells* DNA 
10 repair system in vivo. 

In several methods, the introduction of annealed products into host cells 
selects for heterodupiexes relative to transformed cells comprising homoduplexes. Such 
can be achieved, for example, by providing a first polynucleotide variant as a component 
of a first vector, and a second polynucleotide variant is provided as a component of a 
15 second vector. The first and second vectors are converted to linearized forms in which 
the first and second polynucleotide variants occur at opposite ends. In the incubating 
step, single-stranded forms of the first linearized vector reanneal with each other to form 
linear first vector, single-stranded forms of the second linearized vector reanneal with 
each other to form linear second vector, and single-stranded linearized forms of the first 
20 and second vectors anneal with each to form a circular heteroduplex bearing a nick in 

each strand. Introduction of the products into cells thus selects for cirular heterodupiexes 
relative to the linear first and second vector. Optionally, in the above methods, the first 
and second vectors can be converted to linearized forms by PCR. Alternatively, the first 
and second vectors can be convened to linearized forms by digestion with first and 
25 second restriction enzymes. 

In some methods, polynucleotide variants are provided in double stranded 
form and are converted to single stranded form before the annealing step. Optionally, 
such conversion is bv conducting asymmetric amplification of the first and second 
double stranded polynucleotide variants to amplify a first strand of the first 
30 polynucleotide variant, and a second strand of the second polynucleotide variant. The 
first and second strands anneal in the incubating step to form a heteroduplex. 

in some methods, a population of polynucleotides comprising first and 
second poivnucieotides is provided in double stranded form, and the method further 
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comprises incorporating the first and second polynucleotides as components of first and 
second vectors, whereby the first and second polynucleotides occupy opposite ends of 
the first and second vectors. In the incubating step single-stranded forms of the first 
linearized vector reanneal with each other to form linear first vector, single-stranded 
forms of the second linearized vector reanneal with each other to form linear second 
vector, and single-stranded linearized forms of the first and second vectors anneal with 
each to form a circular heteroduplex bearing a nick in each strand. In the introducing 
step selects for transformed cells comprises the circular heteroduplexes relative to the 
linear first and second vector. 

In some methods, the first and second polynucleotides are obtained from 
chromosomal DNA. In some methods, the polynucleotide variants encode variants of a 
polypeptide. In some methods, the population of polynucleotide variants comprises at 
least 20 variants. In some methods, the population of polynucleotide variants are at least 
1 0 kb in length. 

In some methods, the polynucleotide variants comprises natural variants. 
In other methods, the polynucleotide variants comprise variants generated by mutagenic 
PCR or cassette mutagenesis. In some methods, the host cells into which heteroduplexes 
are introduced are bacterial cells. In some methods, the population of variant 
polynucleotide variants comprises at least 5 polynucleotides having at least 90% sequence 
) identity with one another. 

Some methods further comprise a step of at least partially demethylating 
variant polynucleotides. Demethylation can be achieved by PCR amplification or by 
passaging variants through methylation-deficient host cells. 

Some methods include a further step of sealing one or more nicks in 
5 heteroduplex molecules before exposing the heteroduplexes to a DN A repair system. 
Nicks can be sealed by treatment with DNA ligase. 

Some methods further comprise a step of isolating a screened recombinant 
polynucleotide anant. In some methods, the polynucleotide variant is screened to 
produce a recombinant protein or a secondary metabolite whose production is catalyzed 
0 thereby. 

In some methods, the recombinant protein or secondary 1 metabolite is 
formulated with a carrier to form a pharmaceutical composition. 
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In some methods, the polynucleotide variants encode enzymes selected 
from the group consisting of proteases, lipases, amylases, cutmases. cellulases. amylases, 
oxidases, peroxidases and phytases. In other methods, the polynucleotide variants encode 
a polypeptide selected from the group consisting of insulin. ACTH. glucagon. 

5 somatostatin, somatotropin, thymosin, parathyroid hormone, pigmentary hormones, 

somatomedin, erthropoietin. luteinizing hormone, chorionic gonadotropin, hyperthalmic 
releasing factors, antidiuretic hormones, thyroid stimulating hormone, relaxin. interferon, 
thrombopoietic (TPO), and prolactin. 

In some methods, each polynucleotide in the population of variant 

10 polynucleotides encodes a plurality of enzymes forming a metabolic pathway. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 illustrates the process of heteroduplex formation using 
polymerase chain reaction (PCR) with one set of pnmers for each different sequence to 
1 5 amplify the target sequence and vector. 

Figure 2 illustrates the process of heteroduplex formation using restriction 
enzymes to linearize the target sequences and vector. 

Figure 3 illustrates a process of heteroduplex formation using asymmetric 
or single primer polymerase chain reaction (PCR) with one set of primers for each 
20 different sequence to amplify the target sequence and vector. 

Figure 4 illustrates heteroduplex recombination using unique restriction 
enzymes (X and Y) to remove the homoduplexes. 

Figure 5 shows the amino acid sequences of the FlaA from R. lupini (SEQ 
ID NO: 1) and R. meliloti (SEQ ID NO:2). 
25 Figures 6A and 6B show the locations of the unique restriction sites 

utilized to linearize pRL20 and pRM40. 

Figures 7 A. B. C and D show the DNA sequences of four mosaic J1cl4 
eenes created by in vitro heteroduplex formation followed by in vivo repair ((a) is SEQ 
ID NO;3. (b) is SEQ ID NO:4. (ciis SEQ ID NO:5 and id) is SEQ ID NO:6). 
30 Figure 8 illustrates how the heteroduplex repair process created mosaic 

flaA genes containing sequence information from both parent genes. 
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Ficure 9 shows the physical maps of Actinoplanes utahensis ECB 
deacylase mutants with enhanced specific activity ((a) is pM7-2 for Mutant 7-2, and (b) is 
pM16 for Mutant 16). 

Figure 10 illustrates the process used for Example 2 to recombine 
5 mutations in Mutant 7-2 and Mutant 16 to yield ECB deacylase recombinant with more 
enhanced specific activity. 

Figure 1 1 shows specific activities of wild-type ECB deacylase and 
improved mutants Mutant 7-2, Mutant 16 and recombined Mutant 15. 

Figure 12 shows positions of DNA base changes and amino acid 
1 0 substitutions in recombined ECB deacylase Mutant 1 5 with respect to parental sequences 
of Mutant 7-2 and Mutant 16. 

Figures 13 A, B, C, D and E show the DNA sequence of A.wahensis ECB 
deacylase gene mutant M-15 genes created by in vitro heteroduplex formation followed 
by in vivo repair (SEQ ID NO:7). 
15 Figure 14 illustrates the process used for Example 3 to recombine 

mutations in RC1 and RC2 to yield thermostable subtilisin E. 

Figure 1 5 illustrates the sequences of RC1 and RC2 and the ten clones 
picked randomly from the transformants of the reaction products of duplex formation as 
described in Example 3. The x's correspond to base positions that differ between RC1 and 
20 RC2. The mutation at 995 corresponds to amino acid substitution at 181, while that at 
1 107 corresponds to an amino acid substitution at 21 8 in the subtilisin protein sequence. 

Figure 16 shows the results of screening 400 clones from the library 
created by heteroduplex formation and repair for initial activity (AO and residual activity 
(A r ). The ratio A,/ A r was used to estimate the enzymes' thermostability. Data from 
25 active variants are sorted and plotted in descending order. Approximately 12.9% of the 
clones exhibit a phcnotype corresponding to the double mutant containing both the 
N181D and the N218S mutations. 



DEFINITIONS 

30 Screening is. in general, a two-step process in which one first physically 

separates the cells and then determines which cells do and do not possess a desired 
property. Selection is a form of screening in which identification and physical separation 
are achieved simultaneously by expression of a selection marker, which, in some genetic 
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circumstances, allows ceils expressing the marker to survive while other ceils die t or vice 
versai. Exemplary screening memoers include luciferase, igalactosidase and green 
fluorescent protein. Selection markers include drug and toxin resistance genes. Although 
spontaneous selection can and does occur in the course of natural evolution, in the present 
5 methods selection is performed by man. 

An exogenous DNA segment is one foreign (or heterologous) to the ceil or 
homologous to the cell but in a position within the host cell nucleic acid in which the 
clement is not ordinarily found. Exogenous DNA segments are expressed to yield 
exogenous polypeptides. 
10 The term gene is used broadly to refer to any segment of DNA associated 

with a biological function. Thus, genes include coding sequences and/or the regulatory 
sequences required for their expression. Genes also include nonexpressed DNA segments 
that, for example, form recognition sequences for other proteins. 

The term "wild-type" means that the nucleic acid fragment does not 
1 5 compose any mutations. A "wild-type" protein means that the protein will be active at a 
level of activity found in nature and typically will comprise the amino acid sequence 
found in nature. In an aspect, the term "wild type" or "parental sequence" can indicate a 
starting or reference sequence prior to a manipulation of the invention. 

"Substantially pure" means an object species is the predominant species 
20 present ( i.e., on a molar basis it is more abundant than any other individual 

macromolecular species in the composition), and preferably a substantially purified 
fraction is a composition wherein the object species comprises at least about 50 percent 
(on a molar basis) of all macromolecular species present. Generally, a substantially pure 
composition will comprise more than about 80 to 90 percent of all macromolecular 
25 species present in the composition. Most preferably, the object species is purified to 
essential homogeneity (contaminant species cannot be detected in the composition by 
conventional detection methods) wherein the composition consists essentially of a single 
macromolecular species. Solvent species, small molecules ( <500 Daltons). and elemental 
ion species are not considered macromolecular species. 
30 Percentage sequence identity is calculated by comparing two optimally 

aliened sequences over the window of comparison, determining the number of positions 
at which the identical nucleic acid base occurs in both sequences to yield the number of 
matched positions, dividing the numoer of matched positions by the iota; number o: 



WO 99/29902 



PCT/US98/25698 



- 8 - 

positions in the window of comparison. Optimal alignment of sequences for aligning a 
comparison window can be conducted by computerized implementations of algorithms 
GAP, BESTFIT. FASTA. and TFASTA in the Wisconsin Genetics Software Package 
Release 7.0, Genetics Computer Group, 575 Science Dr.. Madison. WI. 

The term naturally-occurring is used to describe an object that can be 
found in nature as distinct from being artificially produced by man. For example, a 
polypeptide or polynucleotide sequence that is present in an organism (including viruses) 
that can be isolated from a source in nature and which has not been intentionally modified 
by man in the laboratory is naturally-occurring. Generally, the term naturally-occurring 
refers to an object as present in a non-pathological (undiseased) individual, such as would 
be typical for the species. 

A nucleic acid is operably linked when it is placed into a functional 
relationship with another nucleic acid sequence. For instance, a promoter or enhancer is 
operably linked to a coding sequence if it increases the transcription of the coding 
sequence. Operably linked means that the DNA sequences being linked are typically 
contiguous and. where necessary to join two protein coding regions, contiguous and in 
reading frame. However, since enhancers generally function when separated from the 
promoter by several kilobases and intronic sequences may be of variable lengths, some 
polynucleotide elements may be operably linked but not contiguous. 

A specific binding affinity between, for example, a ligand and a receptor, 
means a binding affinity of at least 1 x 10 6 NT 1 . 

The term "cognate" as used herein refers to a gene sequence that is 
evolutionarily and functionally related between species. For example but not limitation, 
in the human genome, the human CD4 gene is the cognate gene to the mouse CD4 gene, 
5 since the sequences and structures of these two genes indicate that they are highly 
homologous and both genes encode a protein which functions in signaling T cell 
activation through MHC class II-restricted antigen recognition. 

The term "heteroduplex" refers to hybrid DNA generated by base pairing 
between complementary single strands derived from the different parental duplex 
0 molecules, whereas the term ^homoduplex" refers to double-stranded DNA generated b 
base pairing between complementary single strands derived from the same parental 
duplex molecules. 
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The term '"nick" in duplex DNA refers to the absence of a phosphochester 
bond between two adjacent nucleotides on one strand. The term "cap" in duplex DNA 
refers to an absence of one or more nucleotides in one strand of the duplex. The term 
Hoop" in duplex DNA refers to one or more unpaired nucleotides in one strand. 
5 A mutant or variant sequence is a sequence showing substantial variation 

from a wild type or reference sequence that differs from the wild type or reference 
sequence at one or more positions. 

DETAILED DESCRIPTION 

10 

I. General 

The invention provides methods of evolving a polynucleotide toward 
acquisition of a desired property. The substrates for the method are a population of at 
least two polynucleotide variant sequences that contain regions of similarity with each 
1 5 other but, which also have point(s) or regions of divergence. The substrates are annealed 
in vitro at the regions of similarity. .Annealing can regenerate initial substrates or can 
form heteroduplexes. in which component strands originate from different parents. The 
products of annealing are exposed to enzymes of a DNA repair, and optionally a 
replication system, that repairs unmatched pairings. Exposure can be in vivo as when 
20 annealed products are transformed into host cells and exposed to the hosts DNA repair 
system. Alternatively, exposure can be in vitro, as when annealed products are exposed 
to cellular extracts containing functional DNA repair systems. Exposure of 
heteroduplexes to a DNA repair system results in DNA repair at bulges in the 
heteroduplexes due to DNA mismatching. The repair process differs from homologous 
25 recombination in promoting nonreciprocal exchange of diversity between strands. The 
DNA repair process is typically effected on both component strands of a heteroduplex 
molecule and at any particular mismatch is typically random as to which strand is 
repaired. The resulting population can thus contain recombinant polynucleotides 
encompassing an essentially random reassortment of points of divergence between 
30 parental strands. The population of recombinant polynucleotides is then screened for 
acquisition of a desired property. The property can be a property of the polynucleotide 
ncr sc. such as capacity of a DNA molecule to bind to a protein or can be a property of an 
expression product thereof, such as mRNA or a protein. 
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II. Substrates For Shuffling 

The substrates for shuffling are variants of a reference polynucleotide that 
show some region(s) of similarity with the reference and other region(s) or pomt(s) of 
divergence. Regions of similarity should be sufficient to support annealing of 
polynucleotides such that stable heteroduplexes can be formed. Variants forms often 
show substantial sequence identity with each other (e.g., at least 50%. 75%, 90% or 99%). 
There should be at least sufficient diversity between substrates that recombination can 
generate more diverse products than there are starting materials. Thus, there must be at 
0 least two substrates differing in at least two positions. The degree of diversity depends on 
the length of the substrate being recombined and the extent of the functional change to be 
evolved. Diversity at between 0.1-25% of positions is typical. Recombination of 
mutations from very closely related genes or even whole sections of sequences from more 
distantly related genes or sets of genes can enhance the rate of evolution and the 
5 acquisition of desirable new properties. Recombination to create chimeric or mosaic 
genes can be useful in order to combine desirable features of two or more parents into a 
single gene or set of genes, or to create novel functional features not found in the parents. 
The number of different substrates to be combined can vary widely in size from two to 
10. 100. 1000, to more than 10 5 , 10 7, or 10 9 members. 
0 The initial small population of the specific nucleic acid sequences having 

mutations may be created by a number of different methods. Mutations may be created 
by error-prone PCR. Error-prone PCR uses low-fidelity polymerization conditions to 
introduce a low level of point mutations randomly over a long sequence. Alternatively, 
mutations can be introduced into the template polynucleotide by oligonucleotide-directed 
15 mutagenesis. In oligonucleotide-directed mutagenesis, a short sequence of the 

polynucleotide is removed from the polynucleotide using restriction enzyme digestion 
and is replaced with a synthetic polynucleotide in which various bases have been altered 
from the original sequence. The polynucleotide sequence can also be altered by chemical 
mutagenesis. Chemical mutagens include, for example, sodium bisulfite, nitrous acid. 
30 hydroxylamme. hydrazine or formic acid. Other agents which are analogues of 

nucleotide precursors include nitrosoguanidme. 5-bromouracil. 2-aminopunne. or 
acridine. Generally, these agents are added to the PCR reaction in place of the nucleotide 
precursor therebv mutating the sequence. Intercalating agents such as proflavine. 
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acrifiavine. qumacrine and the like can also be used. Random mutagenesis of the 
polynucleotide sequence can also be achieved by irradiation with X-rays or ultraviolet 
light. Generally, piasmid DNA or DNA fragments so mutagenized are introduced into E. 
coii and propagated as a pool or library of mutant plasmids. 
5 Alternatively the small mixed population of specific nucleic acids can be 

found in nature in the form of different alleles of the same gene or the same gene from 
different related species ( i.e., cognate genes). Alternatively, substrates can be related but 
nonallelic genes, such as the immunoglobulin genes. Diversity can also be the result of 
previous recombination or shuffling. Diversity can also result from resynthesizing genes 
0 encoding natural proteins with alternative codon usage. 

The starting substrates encode variant forms of sequences to be evolved. 
In some methods, the substrates encode variant forms of a protein for which evolution of 
a new or modified property is desired. In other methods, the substrates can encode 
variant forms of a plurality of genes constituting a multigene pathway. In such methods. 
5 variation can occur in one or any number of the component genes. In other methods, 
substrates can contain variants segments to be evolved as DNA or RNA binding 
sequences. In methods, in which starting substrates containing coding sequences, any 
essential regulatory sequences, such as a promoter and polyadenylation sequence, 
required for expression may also be present as a component of the substrate. 
20 Alternatively, such regulatory sequences can be provided as components of vectors used 
for cloning the substrates. 

The starting substrates can vary in length from about 50, 250, 1000, 
10.000, 100,000, 10 6 or more bases. The starting substrates can be provided in double- or 
single-stranded form. The starting substrates can be DNA or RNA and analogs thereof. 
25 If DNA. the starting substrates can be genomic or cDNA. If the substrates are RNA. the 
substrates are typically reverse-transcribed to cDNA before heteroduplex formation. 
Substrates can be provided as cloned fragments, chemically synthesized fragments or 
PCR amplification products. Substrates can derive from chromosomal, piasmid or viral 
sources. In some methods, substrates are provided in concatemenc form. 

30 

HI. Procedures for Generating Heterodupiexes 

Heteroduplexes are generated from double stranded DNA substrates, by 
denaturing the DNA substrates and incubating under annealing conditions. Hybridization 
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conditions for heteroduplex formation are sequence-dependent and are different in 
different circumstances. Longer sequences hybridize specifically at higher temperatures. 
Generally, hybridization conditions are selected to be about 25°C lower than the thermal 
melting point (Tm) for the specific sequence at a defined ionic strength and pH. The Tm 
is the temperature (under defined ionic strength, pH, and nucleic acid concentration) at 
which 50% of the probes complementary to the target sequence hybridize to the target 
sequence at equilibrium. 

Exemplary conditions for denaturation and renaturation of double stranded 
substrates are as follows. Equimolar concentrations (~ 1 .0 - 5.0 nM) of the substrates are 
mixed in 1 x SSPE buffer ( 1 80 mM NaCl. 1 .0 mM EDTA, 1 0 mM NaH 2 P0 4 , pH 7.4) 
After heating at 96°C for 10 minutes, the reaction mixture is immediately cooled at 0°C 
for 5 minutes; The mixture is then incubated at 68°C for 2-6 hr. Denaturation and 
reannealing can also be carried out by the addition and removal of a denaturant such as 
NaOH. The process is the same for single stranded DNA substrates, except that the 
denaturing step may be omitted for short sequences. 

By appropriate design of substrates for heteroduplex formation, it is 
possible to achieve selection for heteroduplexes relative to reformed parental 
homoduplexes. Homoduplexes merely reconstruct parental substrates and effectively 
dilute recombinant products in subsequent screening steps. In general, selection is 
) achieved by designing substrates such that heteroduplexes are formed in open-circles, 

whereas homoduplexes are formed as linear molecules. A subsequent transformation step 
results in substantial enrichment (e.g., 1 00-fold) for the circular heteroduplexes. 

Figure 1 shows a method in which two substrate sequences in separate 
vectors are PCR-amplified using two different sets of primers (PI, P2 and P3. P4). 
5 Typically, first and second substrates are inserted into separate copies of the same vector. 
The two different pairs of primers initiate amplification at different points on the two 
vectors. Fig. 1 shows an arrangement in which the P1/P2 primer pairs initiates 
amplification at one of the two boundaries of the vector with the substrate and the P1/P2 
primer pair initiates replication at the other boundary in a second vector. The two primers 
0 in each primer pair prime amplification in opposite directions around a circular plasmid. 
The amplification products generated by this amplification are double-stranded linearized 
vector molecules in which the first and second substrates occur at opposite ends of the 
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vector. The amplification products are mixed, denaturec and annealed. Mixing and 
denaturation can be performed in either order. Reannealing generates two linear 
homodupiexes, and an open circular heteroduplex containing one nick in each strand, at 
the initiation point of PCR amplification. Introduction of the amplification products into 
host cells selects for the heteroduplexes relative to the homodupiexes because the former 
transform much more efficiently than the latter. 

It is not essential in the above scheme that amplification is initiated at the 
interface between substrate and the rest of the vector. Rather, amplification can be 
initiated at any points on two vectors bearing substrates provided that the amplification is 
initiated at different points between the vectors. In the general case, such amplification 
generates two linearized vectors in which the first and second substrates respectively 
occupy different positions relative to the remainder of the vector. Denaturation and 
reanneaiing generator heteroduplexes similar to that shown in Fig. 1. except that the nicks 
occur within the vector component rather than at the interface between plasmid and 
substrate. Initiation of amplification outside the substrate component of a vector has the 
advantage that it is not necessary to design primers specific for the substrate borne by the 
vector. 

Although Fig. 1 is exemplified for two substrates, the above scheme can be 
extended to any number of substrates. For example, an initial population of vector 
bearing substrates can be divided into two pools. One pool is PCR-amplified from one 
set of primers, and the other pool from another. The amplification products are denatured 
and annealed as before. Heteroduplexes can form containing one strand from any 
substrate in the first pool and one strand from any substrate in the second pool. 
Alternatively, three or more substrates cloned into multiple copies of a vector can be 
5 subjected to amplification with amplification in each vector starting at a different point. 
For each substrate, this process generates amplification products varying in how flanking 
vector DNA is divided on the two sides of the substrate. For example, one amplification 
product might have most of the vector on one side of the substrate, another amplification 
product might have most of the vector on the other side of the substrate, and a further 
0 amplification product might have an equal division of vector sequence flanking the 
substrate. In the subsequent annealing step, a strand of substrate can form a circular 
heteroduplex with a strand of any other substrate, but strands of the same substrate can 
oniv reanneai with each other to form a linear homodupiex. In a still further variation. 
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multiple substrates can be performed by performing multiple iterations of the scheme in 
Fig. 1 . After the first iteration, recombinant polynucleotides in a vector, undergo 
heteroduplex formation with a third substrate incorporated into a further copy of the 
vector. The vector bearing the recombinant polynucleotides and the vector bearing the 
5 third substrate are separately PCR amplified from different primer pairs. The 

amplification products are then denatured and annealed. The process can be repeated 
further times to allow recombination with further substrates. 

An alternative scheme for heteroduplex formation is shown in Fig. 2. 
Here, first and second substrates are incorporated into separate copies of a vector. The 
10 two copies are then respectively digested with different restriction enzymes. Fig. 2 shows 
an arrangement in which, the restriction enzymes cut at opposite boundaries between 
substrates and vector, but all that is necessary is to use two different restriction enzymes 
that cut at different places. Digestion generates linearized first and second vector bearing 
first and second substrates, the first and second substrates occupying different positions 
1 5 relative to the remaining vector sequences. Denaturation and reannealing generates open 
circular heteroduplexes and linear homoduplexes. The scheme can be extended to 
recombination between more than two substrates using analogous strategies to those 
described with respect to Fig. 1 . In one variation, two pools of substrates are formed, and 
each is separately cloned into vector. The two pools are then cute with different enzymes, 
20 and annealing proceeds as for two substrates. In another variation, three or more 

substrates can be cloned into three or more copies of vector, and the three or more result 
molecules cut with three or more enzymes, cutting at three or more sites. This generates 
three different linearized vector forms differing in the division of vector sequences 
flanking the substrate moiety in the vectors. Alternatively, any number of substrates can 
25 be recombined pairwise in an iterative fashion with products of one round of 
recombination annealing with a fresh substrate in each round. 

In a further variation, heterodupiexes can be formed from substrates 
molecules in vector-free form, and the heteroduplexes subsequently cloned into vectors. 
Such can be achieved by asymmetric amplification of first and second substrates as 
30 shown in Fig. 3. Asymmetric or single primer PCR amplifies only one strand of a duplex. 
By appropriate selection of primers, opposite strands can be amplified from two different 
substrates. On reannealing amplification products, heteroduplexes are formed from 
opposite strands of the two substrates. Because only one strand is amplified from each 
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substrate, reanneaiing does not reform homodupiexes (other than for small quantities of 
unamplified substrate ). The process can be extended to allow recombination of any 
number of substrates using analogous strategies to those described with respect to Fig. 1. 
For example, substrates can be divided into two poois. and each pool subject to the same 
5 asymmetric amplification, such that amplification products of one pool can only anneal 
with amplification products of the other pool . and not with each other. Alternatively, 
shuffling can proceed pairwise in an iterative manner, in which recombinants formed 
from heteroduplexes of first and second substrates, are subsequently subjected to 
heteroduplex formation with a third substrate. Point mutations can also be introduced at a 
10 desired level during PCR amplification. 

Fig. 4 shows another approach of selecting for heteroduplexes relative to 
homodupiexes. First and second substrates are isolated by PCR amplification from 
separate vectors. The substrates are denatured and allowed to anneal forming both 
heteroduplexes and reconstructed homodupiexes. The products of annealing are digested 
1 5 with restriction enzymes X and Y. X has a site in the first substrate but not the second 
substrate, and vice versa for Y. Enzyme X cuts reconstructed homoduplex from the first 
substrate and enzyme Y cuts reconstructed homoduplex from the second substrate. 
Neither enzyme cuts heteroduplexes. Heteroduplexes can effectively be separated from 
restriction fragments of homodupiexes by further cleavage with enzymes A and B having 
20 sites proximate to the ends of both the first and second substrates, and ligation of the 

products into vector having cohesive ends compatible with ends resulting from digestion 
with A and B. Only heteroduplexes cut with A and B can ligate with the vector. 
Alternatively, heteroduplexes can be separated from restriction fragments of 
homodupiexes by size selection on gels. The above process can be generalized to N 
25 substrates by cleaving the mixture of heteroduplexes and homodupiexes with N enzymes, 
each one of which cuts a different substrate and no other substrate. Heteroduplexes can 
be formed by directional cloning. Two substrates for heteroduplex formation can be 
obtained by PCR amplification of chromosomal DNA and joined to opposite ends of a 
linear vector. Directional cloning can be achieved by digesting the vector with two 
30 different enzymes, and digesting or adapting first and second substrates to be respectively 
compatible with cohesive ends of only of the two enzymes used to cut the vector. The 
first and second substrates can thus be Heated at opposite ends of a linearized vector 
fragment. This scheme can be extended to any number of substrates by using principles 
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analogous to those described for Fig. 1. For example, substrates can be divided into two 
pools before ligation to the vector. Alternatively, recombinant products formed by 
heteroduplex formation of first and second substrates, can subsequently undergo 
heteroduplex formation with a third substrate. 

5 

IV. Vectors and Transformation 

In general, substrates are incorporated into vectors either before or after 
the heteroduplex formation step. A variety of cloning vectors typically used in genetic 
engineering are suitable. 
1 o The vectors containing the DN A segments of interest can be transferred 

into the host cell by standard methods, depending on the type of cellular host. For 
example, calcium chloride transformation is commonly utilized for prokaryotic cells, 
whereas calcium phosphate treatment. Lipofection, or electroporation may be used for 
other cellular hosts. Other methods used to transform mammalian cells include the use of 
1 5 Polybrene, protoplast fusion, liposomes, electroporation, and microinjection, and 

biolisitics (see, generally, Sambrook et aU supra). Viral vectors can also be packaged in 
vitro and introduced by infection. The choice of vector depends on the host cells. In 
general, a suitable vector has an origin of replication recognized in the desired host cell, a 
selection maker capable of being expressed in the intended host cells and/or regulatory 
20 sequences to support expression of genes within substrates being shuffled. 

V. Types of Host Cells 

In general any type of cells supporting DNA repair and replication of 
heteroduplexes introduced into the cells can be used. Cells of particular interest are the 

25 standard cell types commonly used in genetic engineering, such as bacteria, particularly, 
E. coli (16, 1 7). Suitable E. coli strains include E. coli mutS* mutL. dam', and/or recA*. 
E.coli XL-10-Gold {[Tef ' Aimer A) 1 83 A(mcrCB-hsdSMR-mrr) 1 '3 encLAl supE44 thi-1 
rccAl gy>rA96 relAl lac Hte ] [F % pro A B lacPZAM15 TnJO (Tef) Amy Cam J). E. coli 
ES1301 muiS [Genotype. lacZ53. mulSlOl : :Tn5, thyA36, rha-5. meiBL deoC. INirrnD- 

30 rrnE)} (20. 24, 28-42). Preferred E. coli strains are E.coli SCSI 1 0 [Genotype: rpsi 
(Sir), thr. leu, enda, thi-L lacy, gaik, gait, ara tona. tsx, dam, dcm. supE44. Ailac- 
proAB). [F, traD36, proA~ B~ lacVZASil 5\ which have normal cellular mismatch repair 
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systems (\7). This strain type repairs mismatches and unmatches in the heterociupiex 
with linle strand-specific preference. Further, because this strain is aam and dcrri . 
piasmid isolated from the strain is unmethyiated and therefore particularly amenable for 
further rounds of DNA duplex formation' mismatch repair (see below). Other suitable 
5 bacterial cells include gram-negative and gram-positive, such as Baciliuy Pseudomonas. 
and Salmonella. 

Eukaryotic organisms are also able to carry out mismatch repair (43-48). 
Mismatch repair systems in both prokaryotes and eukaryotes are thought to play an 
important role in the maintenance of genetic fidelity during DNA replication. Some of 
1 0 the genes that play important roles in mismatch repair in prokaryotes. particularly mutS 

and mutL. have homologs in eukaryotes. in the outcome of genetic recombinations, and in 
genome stability. Wild-type or mutant 5. ccrevisiae has been shown to carry out 
mismatch repair of heteroduplexes (49-56). as have COS-1 monkey cells (57). Preferred 
strains of yeast are Picchia and Saccharomyces. Mammalian cells have been shown to 
1 5 have the capacity to repair G-T to G-C base pairs by a short-patch mechanism (38. 58- 

63). Mammalian cells (e.g., mouse, hamster, primate, human), both cell lines and primary 
cultures can also be used. Such cells include stem cells, including embryonic stem cells, 
zygotes, fibroblasts, lymphocytes. Chinese hamster ovary (CHO), mouse fibroblasts 
(MIH3T3), kidney, liver, muscle, and skin cells. Other eucaryotic cells of interest include 
20 plant cells, such as maize, rice, wheat, cotton, soybean, sugarcane, tobacco, and 

arabidopsis; fish, algae, fungi (aspergillus. podospora. neurosporai. insect (e.g.. baculo 
lepidoptera) (see, Winnacker. "From Genes to Clones." VCH Publishers. N.Y.. (1987), 
which is incorporated herein by reference). 

In vivo repair occurs in a wide variety of prokaryotic and eukaryotic cells. 
25 Use of mammalian cells is advantage in certain application in which substrates encode 

polypeptides that are expressed only in mammalian cells or which are intended for use in 
mammalian cells. However, bacterial and yeast cells are advantageous for screening 
laree libraries due to the higher transformation frequencies attainable in these strains. 

30 V. In Vitro DNA Repair Svsiems 

As an alternative to introducing annealed products into host ceiis. annealed 
products can be exposed a DNA repair system m vara. The DNA repair system can be 
obtained as extracts from repair-competent E. cuii. yeast or any other ceiis < 64-67V 
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Repair-competent cells are lysed in appropriate buffer and supplemented with 
nucleotides, DNA is incubated in this cell extract and transformed into competent cells 
for replication. 

5 VI. Screening and Selection 

After introduction of annealed products into host cells, the host cells are 
typically cultured to allow repair and replication to occur and optionally, for genes 
encoded by polynucleotides to be expressed. The recombinant polynucleotides can be 
subject to further rounds of recombination using the heteroduplex procedures described 
1 0 above, or other shuffling methods described below. However, whether after one cycle of 
recombination or several, recombinant polynucleotides are subjected to screening or 
selection for a desired property. In some instances, screening or selection in performed in 
the same host cells that are used for DNA repair. In other instances, recombinant 
polynucleotides, their expression products or secondary metabolites produced by the 
1 5 expression products are isolated from such cells and screened in vitro. In other instances, 
recombinant polynucleotides are isolated from the host cells in which recombination 
occurs and are screened or selected in other host cells. For example, in some methods, it 
is advantageous to allow DNA repair to occur in a bacterial host strain, but to screen an 
expression product of recombinant polynucleotides in eucaryotic cells. The recombinant 
20 polynucleotides surviving screening or selection are sometimes useful products in 

themselves. In other instances, such recombinant polynucleotides are subjected to further 
recombination with each other or other substrates. Such recombination can be effected by 
the heteroduplex methods described above or any other shuffling methods. Further 
round(s) of recombination are followed by further rounds of screening or selection on an 
25 iterative basis. Optionally, the stringency of selection can be increased at each round. 

The nature of screening or selection depends on the desired property 
sought to be acquired. Desirable properties of enzymes include high catalytic activity, 
capacity to confer resistance to drugs, high stability, the ability to accept a wider (or 
narrower) range of substrates, or the ability to function in nonnatural environments such 
30 as organic solvents. Other desirable properties of proteins include capacity to bind a 
selected target, secretion capacity, capacity to generate an immune response to a given 
tarcet. lack of immunogenicity and toxicity to pathogenic microorganisms. Desirable 
properties of DNA or RNA polynucleotides sequences include capacity to specifically 
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bind a given protein target, and capacity to regulate expression of operabiy linked coding 
sequences. Some of the above properties, such as drug resistance, can be selected by 
plating cells on the drug. Other properties, such as the influence of a regulatory sequence 

n expression, can be screened by detecting appearance of the expression product of a 
reporter gene linked to the regulator.- sequence. Other properties, such as capacity of an 
expressed protein to be secreted, can be screened by FACS™. using a labelled antibody to 
the protein. Other properties, such as immunoeemciry or lack thereof, can be screened by 
isolating protein from individual cells or pools of cells, and analyzing the protein in vitro 
or in a laboratorv animal. 



VII. Variations 

1. Demethvlation 

Most cell types methylate DNA in some manner, with the pattern of 
methylation differing between cells types. Sites of methylation include 5-methylcytosine 
5 (nrO. N4-methylcytosine (m 4 0 and N 6 -methyladenine (m 6 A\ 5 -hydroxy methylcytosine 
(hm 5 C) and 5-hydroxymethyluracil (WlA In E. coli. methylation is effected by Dam 
and Dcm enzymes. The methylase specified by the dam gene mesylates the N6-position 
of the adenine residue in the sequence GATC. and the methylase specified by the dcm 
gene methylates the C5-position of the internal cytosine residue in the sequence 
10 CCWGG. DNA from plants and mammal is often subject to CG methylation meaning 
that CGorCNG sequences are methylated. Possible effects of methylated on cellular 
repair are discussed by references 18-20. 

In some methods. DNA substrates for heteroduplex formation are at least 
partially demethylated on one or both strands, preferably the latter. Demethvlation of 
25 substrate DNA promotes efficient and random repair of the heteroduplexes. In 

heteroduplexes formed with one strand dam-methylated and one strand unmethylated. 
repair is biased to the unmethylated strand, with the methylated strand serving as the 
template for correction. If neither strand is methylated, mismatch repair occurrs. but 
showes insignificant strand preference (23. 24). 
30 Demethvlation can be performed in a variety of ways. In some methods, 

substrate DNA is demethylated by PCR-ampiification. In some instances. DNA 
demethvlation is accomplished in one of the PCR steps in the heterodupiex formation 
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procedures described above. In other methods, an additional PCR step is performed to 
effect demethylation. In other methods, demethylation is effected by passaging substrate 
DNA through methylation deficient host cells (e.g. an E. coli darridcm strain). In other 
methods, substrate DNA is demethylated in vitro using a demethylating enzymes. 
Demethylated DNA is used for heteroduplex formation using the same procedures 
described above. Heteroduplexes are subsequently introduced into DNA-repair-proficient 
but restriction-enzyme-defective cells to prevent degradation of the unmethylated 
heteroduplexes. 

2. Sealing Nicks 

Several of the methods for heteroduplex formation described above result 
in circular heteroduplexes bearing nicks in each strand. These nicks can be sealed before 
introducing heteroduplexes into host cells. Sealing can be effected by treatment with 
DNA ligase under standard ligating conditions. Ligation forms a phosphodiester bond to 
link two adjacent bases separated by a nick in one strand of double helix of DNA. 
Sealing of nicks increases the frequency of recombination after introduction of 
heteroduplexes into host cells. 

3. Error Prone PCR Attendant To Amplification 

Several of the formats described above include a PCR amplification step. 
Optionally, such a step can be performed under mutagenic conditions to induce additional 
diversity between substrates. 

VIII. Other Shuffling Methods 

The methods of heteroduplex formation described above can be used in 
conjunction with other shuffling methods. For example, one can perform one cycle of 
heteroduplex shuffling, screening or selection, followed by a cycle of shuffling by another 
method, followed by a further cycle of screening or selection. Other shuffling formats are 
described by WO 95/22625; US 5.605.793; US 5.81 1 .238; WO 96/19256; Stemmer. 
Science 270. 1510 (1 995'); Stemmer et al.. Gene, 1 64. 49-53 1 1 995); Stemmer. 
Bio/Technology, 13. 549-553 (1995); Stemmer. Proc. Natl. Acad. Sci. USA 91. 10747- 
10751 (1994); Stemmer. Nature 370. 389-391 (1994); Cramen et al.. Nature Medicine. 
2(1): 1-3, (1996); Cramen et al.. Nature Biotechnology-- 14. 315-319 d 996): WO 9S/4272' 7 : 
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WO 98 41622; WO 98/05764 and WO 98'4272S. WO 98/27230 (each of which is 
incorporated by reference in its entirety for all purposes). 

IX. Protein Anaiogs 

5 Proteins isolated by the methods also serve as lead compounds for the 

development of derivative compounds. The derivative compounds can include chemical 
modifications of amino acids or replace amino acids with chemical structures. The 
anaiogs should have a stabilized electronic configuration and molecular conformation that 
allows key functional groups to be presented in substantially the same way as a lead 
10 protein. In particular, the non-peptic compounds have spatial electronic properties which 
are comparable to the polypeptide binding region, but will typically be much smaller 
molecules than the polypeptides, frequently having a molecular weight below about 2 
CHD and preferably below about 1 CUD. Identification of such non-peptic compounds 
can be performed through several standard methods such as self-consistent field (CSF) 

15 analysis, configuration interaction (CHI) analysis, and normal mode dynamics analysis. 
Computer programs for implementing these techniques are readily available. See Rein et 
al.. Computer-Assisted Modeling of Receptor -Ligand Interactions (Alan Liss, New York. 
1989). 

20 IX. Pharmaceutical Compositions 

Polynucleotides, their expression products, and secondary metabolites 
whose formation is catalyzed by expression products, generated by the above methods are 
optionally formulated as pharmaceutical compositions. Such a composition comprises 
one or more active agents, and a pharmaceuticaliy acceptable earner. A variety of 

25 aqueous carriers can be used. e.g.. water, buffered water, phosphate-buffered saline 

(PBS), 0.4% saline, 0.3% glycine, human albumin solution and the like. These solutions 
are sterile and generally free of particulate matter. The compositions may contain 
pharmaceuticaliy acceptable auxiliary substances as required to approximate physio- 
logical conditions such as pH adjusting and buffering agents, toxicity adjusting agents 

30 and the like, for exampie. sodium acetate, sodium chloride, potassium chloride, calcium 
chloride and sodium is selected primarily based on fluid volumes, viscosities, and so 
forth, in accordance with the particular mode of administration selected. 
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EXAMPLES 

EXAMPLE 1 . Novel Rhizobtum Flag Genes From Recombination Of Rhizohium Lupim 
Flag And Rhizohium Meliloti FlaA 

Bacterial flagella have a helical filament, a proximal hook and a basal 
5 body with the flagellar motor (68). This basic design has been extensively examined in 
£. coli and S. typhimurium and is broadly applicable to many other bacteria as well as 
some archaea. The long helical filaments are polymers assembled from flagellin subunits, 
whose molecular weights range between 20,000 and 65,000, depending on the bacterial 
species (69V Two types of flagellar filaments, named plain and complex, have been 
1 0 distinguished by their electron microscopically determined surface structures (70). Plain 
filaments have a smooth surface with faint helical lines, whereas complex filaments 
exhibit a conspicuous helical pattern of alternating ridges and grooves. These 
characteristics of complex flagellar filaments are considered to be responsible for the 
brittle and (by implication) rigid structure that enables them to propel bacteria efficiently 
1 5 in viscous media (7 1 -73). Whereas flagella with plain filaments can alternate between 

clockwise and counter clockwise rotation (68), all known flagella with complex filaments 
rotate only clockwise with intermittent stops (74). Since this latter navigation pattern is 
found throughout bacteria and archaea, it has been suggested that complex flagella may 
reflect the common background of an ancient, basic motility design (69). 
20 Differing from plain bacterial flagella in the fine structure of their 

filaments dominated by conspicuous helical bands and in their fragility, the filaments are 
also resistant against heat decomposition (72). Schmitt et al. (75) showed that 
bacteriophage 7-7-1 specifically adsorbs to the complex flagella of RAupini H 13-3 and 
requires motility for a productive infection of its host. Though the flagellins from R. 
25 meliloti and R. lupini are quite similar, bacteriophage 7-7-1 does not infect R.meliloti. 
Until now complex flagella have been observed in only three species of soil bacteria: 
Pseudomonas rhodos (73). R.meliloti (76). and RAupini HI 3-3 (70, 72). Cells of RAupini 
HI 3-3 posses 5 to 10 pentrichously inserted complex flagella. which were first isolated 
and analyzed by high resolution electron microscopy and by optical diffraction (70). 
30 Maruyama et al. (77) further found that a higner content of hydrophobic 

amino acid residues in the complex filament may be one of the main reasons for the 
unusual properties of complex flagella. By measuring mass per unit length and obtaining 
three-dimensional reconstruction from electron micrographs. Trachtenberc et al. (73. 78) 
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su£Z2ested that the compiex filaments of R. lupin: are composed of functional aimers. 
Fieure 6 shows the comparison between the deduced ammo acid sequence of the R. lupim 
HI 3-3 FlaA and the deduced amino acid sequence of the R. meliioti FlaA. Perfect 
matches are indicated by vertical lines, and conservative exchanges are indicated by 
5 colons. The overall identity is 56%. The R.lupimflaA and R.meiiloti flaA were subjected 
to in vitro heteroduplex formation followed by in vivo repair in order to create novel 
FlaA molecules and structures. 

A. Methods 

0 pRL20 containing R. lupini HI 3-3 flaA gene and pRM40 containing 

R.meiiloti flaA gene are shown in Figs. 6A and 6B. These plasmids were isolated from £ 
coli SCSI 10 (free from dam- and dcm-type methylation). 

About 3.0 pg. of unmeihylated pRL20 and pRM40 DNA were digested with Bam HI and 
Eco RI. respectively, at 37°C for 1 hour. After agarose gel separation, the linearized 

5 DNA was purified with Wizard PCR Prep kit (Promega, WI, USA). 

Equimolar concentrations (2.5 nM) of the linearized unmethylated pRL20 and pRM40 
were mixed in 1 x SSPE buffer (1 80 mM NaCl, 1 mM EDTA, 10 mM NaH2P04, pH 
7.4). After heating at 96°C for 10 minutes, the reaction mixture was immediately cooled 
at 0°C for 5 minutes. The mixture was incubated at 68°C for 2 hour for heteroduplexes to 

20 form. 

One microliter of the reaction mixture was used to transform 50 ul of E. 
coli ES 1301 mutS. E. coli SCSI 10 and E. coli JM109 competent cells. The 
transformation efficiency with £. coli JM109 competent cells was about seven times 
higher than that of E. coli SCSI 10 and ten times higher than that of E. coli ESI 301 mutS. 
25 although the overall transformation efficiencies were 10-200 times lower than those of 
control transformations with the close, covalent and circular pUC19 plasmid. 

Two clones were selected at random from the £. coli SCSI 10 
transformants and two from E. coli ESI 301 mutS transformants. and plasmid DNA was 
isolated from these four clones for further DNA sequencing analysis. 
30 B. Results 

Ficure 7 shows (a) the sequence of SCS01 (clone-1 from E. coli SCSI 10 
transformant library i. (b> the sequence of SCS02 (clone p2 from E. coli SCSI 10 
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transibirnant library), (c) the sequence of ES01 (clone #1 from E. co//ES1301 
trans formant library)* and (d) the sequence of ES02 (clone #2 from E. coli ESI 301 
transformant library). All four sequences were different from wild-type R. lupiniflaA and 
R. melilotiflaA sequences. Clones SCS02. ES01 and ES02 all contain a complete open- 
reading frame, but SCS01 was truncated. Figure 8 shows that recombination mainly 
occurred in the loop regions (unmatched regions). The flaA mutant library generated 
from R. meliloti flaA and R. lupiniflaA can be transformed into E. coli SCSI 10, ESI 301, 
XLIO-Gold and JM109. and transformants screened for functional FlaA recombinants. 

EXAMPLE 2. Directed Evolution Of ECB Deacvlase For Variants With Enhanced 
Specific Activity 

Streptomvces are among the most important industrial microorganisms due 
to their ability to produce numerous important secondary metabolites (including many 
antibiotics) as well as large amounts of enzymes. The approach described here can be 
used with little modification for directed evolution of native Streptomvces enzymes, 
some or all of the genes in a metabolic pathways, as well as other heterologous enzymes 
expressed in Streptomyces. 

New antifungal agents are critically needed by the large and growing 
numbers of immune-compromised AIDS, organ transplant and cancer chemotherapy 
patients who suffer opportunistic infections. Echinocandin B (ECB), a lipopeptide 
produced by some species of Aspergillus, has been studied extensively as a potential 
antifungal. Various antifungal agents with significantly reduced toxicity have been 
generated by replacing the linoieic acid side chain of A. nidulans echinocandin B with 
different aryl side chains (79-83). The cyclic hexapeptide ECB nucleus precursor for the 
5 chemical acylation is obtained by enzymatic hydrolysis of ECB using Actinoplanes 

utahensis ECB deacvlase. To maximize the conversion of ECB into intact nucleus, this 
reaction is earned out at pH 5.5 with a small amount of miscible organic solvent to 
solubilize the ECB substrate. The product cyclic hexapeptide nucleus is unstable at pH 
above 5.5 during the long incubation required to fully deacylate ECB (84). The pH 
0 optimum of ECB deacvlase. however, is 8.0-8.5 and its activity is reduced at pH 5.5 and 
in the presence of more than 2.5% ethanol (84). To improve production of ECB nucleus 
it is necessary to increase the activity of the ECB deacvlase under these process-relevant 
conditions. 
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Relativeiv iittie is known about EC3 deacyiase. The enzyme is a 
heterodimer whose two subunits are derived by processing of a single precursor protein 
<83V The 19.9 kD a-subunit is separated from the 60.4 kD B-subunit by a 15-amino acid 
spacer peptide that is removed along with a signal peptide and another spacer peptide in 
5 the native organism. The polypeptide is also expressed and processed into functional 

enzyme in Streptomyces lividans. the organism used for iarge-scale conversion of ECB by 
recombinant ECB deacyiase. The three-dimensional structure of the enzyme has not been 
determined, and its sequence shows so little similarity to other possibly related enzymes 
such as penicillin acylase that a structural model reliable enough to guide a rational effort 
0 to engineer the ECB deacyiase will be difficult to build. We therefore decided to use 
directed evolution (85) to improve this important activity. 

Protocols suitable for mutagenic PCR and random-priming recombination 
of the 2.4 kb ECB deacyiase gene (73% G+C) have been described recently (86). Here, 
we further describe the use of heteroduplex recombination to generate new ECB 
5 deacyiase with enhanced specific activity. 

In this case, two Actinopianes uiahensis ECB deacyiase mutants, M7~2 
and Ml 6, which show higher specific activity at pH 5.5 and in the presence of 10% 
MeOH were recombined using technique of the in vitro heteroduplex formation and in 
vivo mismatch repair . 

10 Fieure 1 2 shows the physical maps of plasmids pM7-2 and pM 1 6 which 

contain the genes for the M7-2 and M 16 ECB deacyiase mutants. Mutant M7-2 was 
obtained through mutagenic PCR performed directly on whole Streptomyces lividans cells 
containing wild-type ECB deacyiase gene, expressed from plasmid pSHP 150-2*. 
Streptomyces with pM7-2 show 1.5 times the specific activity of cells expressing the 

25 wild-type ECB deacyiase (86). Clone pM16 was obtained using the random-pnming 

recombination technique as described (86. 87). It shows 2.4 times specific activity of the 
wild-type ECB deacyiase clone. 

A. Methods: 

30 M7-2 and M16 plasmid DNA <pM7-2 and pM16) (Fig. 9) were purified 

from E. coil SCSI 10 (in separate reactions). About 5.0 ug of unmethyiated M7-2 and 
Ml 6 DNA were digested with A no 1 and Psh AL respectively, at 37°C for 1 hour (Fig. 
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] 0). After agarose gel separation, the linearized DNA was purified using a Wizard PCR 
Prep Kit (Promega, WI. USA). 

Equimolar concentrations (2.0 nM ) of the linearized unmethylated pM7-2 and pM16 
DNA were mixed in 1 x SSPE buffer ( 1 x SSPE: 1 80 mM NaCl, 1 .0 mM EDTA, 1 0 mM 
5 NaEbPOa, pH 7.4). After heating at 96 °C for 10 minutes, the reaction mixture is 
immediately cooled at 0 °C for 5 minutes. The mixture was incubated at 68 °C for 3 
hours to promote formation of heteroduplexes. 

One microliter of the reaction mixture was used to transform 50 |_il of 
Ecoli ESI 301 mutS. SCSI 10 and JM109 competent cells. All transformants from E. coli 
10 ES1301 mutS were pooled and E. coli SCSI 10 were pooled. A plasmid pool was isolated 
from each pooled library, and this pool w r as used to transform S. lividans TK23 
protoplasts to form a mutant library for deacylase activity screening. 
Transformants from the S. lividans TK23 libraries were screened for ECB deacylase 
activity with an in situ plate assay. Transformed protoplasts were allowed to regenerate 
1 5 on R2YE agar plates for 24 hr at 30°C and to develop in the presence of thiostrepton for 
48 hours. When the colonies grew to the proper size, 6 ml of 0.7% agarose solution 
containing 0.5 mg/ml ECB in 0, 1 M sodium acetate buffer (pH 5.5) was poured on ton of 
each R2YE-agar plate and allowed to develop for 18-24 hr at 30°C. Colonies surrounded 
by a clearing zone larger than that of a control colony containing wild-type plasmid 
20 pSHP150-2*, were selected for further characterization. 

Selected transformants were inoculated into 20 ml medium containing 
thiostrepton and grown aerobically at 30°C for 48 hours, at which point they were 
analyzed for ECB deacylase activity using HPLC. 100 |J.l of whole broth was used for a 
reaction at 30 °C for 30 minutes in 0.1 M NaAc buffer (pH 5.5) containing 10% (v/v) 
25 MeOH and 200 |ig/ml of ECB substrate. The reactions were stopped by adding 2.5 

volumes of methanol, and 20 fj.1 of each sample were analyzed by HPLC on a 100 x 4.6 
mm polyhydroxyethyl aspartamide column (PolyLC Inc.. Columbia, MD, USA) at room 
temperature using a linear acetonurile gradient starting with 50:50 of A:B (A = 93% 
acetonitrile. 0.1% phosphoric acid: B = 70% acetonitrile. 0.1% phosphoric acid) and 
30 ending with 30:70 of A:B in 22 mm at a flow rate of 2.2 ml/min. The areas of the ECB 
and ECB nucleus peaks were calculated and subtracted from the areas of the 
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corresponding peaks from a sample culture of 5. lividans containing pIJ702* in order to 
estimate the ECB deacyiase activity. 

2.0 ml pre-cuitures of positive mutants were used to inoculate 50-ml 
medium and allowed to grow at 30°C for 96 hr. The supernatants were further 
concentrated to 1/30 their original volume using an Amicon filtration unit (Beverly. MA, 
USA) with molecular weight cutoff of 10 kD. The resulting enzyme samples were 
diluted with an equal volume of 50 mM KH 2 P04 (pH 6.0) buffer and were applied to Hi- 
Trap ion exchange column (Pharmacia Biotech, Piscataway, NL USA) . The binding 
buffer was 50 mM KH : P0 4 (pH 6.0). and the elution buffer was 50 mM KH 2 P0 4 (pH 6.0) 
containing 1 .0 M NaCl. A linear gradient from 0 to 1.0 M NaCl was applied in 8 column 
volumes with a flow rate of 2.7 ml/min. The ECB deacyiase fraction eluting at 0.3 M 
NaCl was concentrated and the buffer was exchanged for 50 mM KH 2 P0 4 (pH 6.0) using 
Centncon-10 units. Enzyme purity was verified by SDS-PAGE using Coomassie Blue 
stain, and the concentration was determined using the Bio-Rad Protein Assay Reagent 

(Hercules, CA. USA). 

A modified HPLC assay was used to determine the activities of the ECB 
deacyiase mutants on ECB substrate (84). Four >ig of each purified ECB deacyiase 
mutant was used for activity assay reaction at 30°C for 30 minutes in 0.1 M NaAc buffer 
(pH 5.5) containing 10% (v/v) MeOH and different concentrations of ECB substrate. 
Assays were performed in duplicate. The reactions were stopped by adding 2.5 volumes 
of methanol, and the HPLC assays were earned out as described above. The absorbance 
values were recorded, and the initial rates were calculated by least-squares regression of 
the time progress curves from which the Km and the kcat were calculated. 

Activities as a function of pH were measured for the purified ECB 
deacylases at 30°C at different pH values: 5, 5.5 and 6 (0.1 M acetate buffer); 7. 7.5, 8 
and 8.5 (0.1 M phosphate buffer); 9 and 10 (0.1 M carbonate buffer ) using the HPLC 
assay. Stabilities of purified ECB deacylases were were determined at 30°C in 0.1 M 
NaAc buffer (pH 5.5) containing 10% methanol. Samples were withdrawn at different 
time intervals, and the residual activity was measured in the same buffer with the HPLC 
assav described above. 
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B. Results 

Fig. 1 1 shows that after one round of applying this heteroduplex repair 
technique on the mutant M7-2 and Ml 6 genes, one mutant (Ml 5) from about 500 original 
transformants was found to possess 3.1 times the specific activity of wild-type. 
Wild type and evolved Ml 5 ECB deacylases were purified and their kinetic parameters 
for deacylation of ECB were determined by HPLC. The evolved deacylases Ml 5 has an 
increased catalytic rate constant, k^ai by 205%. The catalytic efficiency (kcat/K m ) of M20 
is enhanced by a factor of 2.9 over the wild-type enzyme. 

Initial rates of deacylation with the wild type and Ml 5 at different pH 
values from 5 to 10 were determined at 200 ng/ml of ECB. The recombined Ml 5 is 
more active than wild type at pH 5-8. Although the pH dependence of the enzyme 
activity in this assay is not strong, there is a definite shift of 1 .0-1 .5 units in the optimum 
to lower pH, as compared to wild type. 

The time courses of deactivation of the purified ECB deacylase mutant 
Ml 5 was measured in 0.1 M NaAc (pH 5.5) at 30°C. No significant difference in 
stability was observed between wild type and mutant Ml 5. 

The DNA mutations with respect to the wild type ECB deacylase sequence 
and the positions of the amino acid substitutions in the evolved variants M7-2, Ml 6 and 
Ml 5 are summarized in Figure 12. 
0 The heteroduplex recombination technique can recombine parent 

sequences to create novel progeny. Recombination of the M7-2 and M16 genes yielded 
M15„ whose activity is higher than any of its parents (Fid. 13). Of the six base 
substitutions in M15, five (at positions a50. al71. 357, pi29 and p340) were inherited 
from M7-2, and the other one (P30) came from Ml 6. 
5 This approach provides an alternative to existing methods of DNA 

recombination and is particularly useful in recombining large genes or entire operons. 
This method can be used to create recombinant proteins to improve their properties or to 
study structure-function relationship. 



W O 99 7990: 



PCT.TS98 '25698 



- 29 - 

EXAMPLE 3 Novel Thermostable Bacillus Suhritis Subtiiisin E Variants 

This example demonstrates the use in vitro heterodupiex ibrmauon 
followed bv in vivo repair tor combining sequence information from two different 
sequences in order to improve the thermostability of Bacillus subtilis subtiiisin E 
5 Genes RC1 and RC2 encode thermostable B. sublilis subtiiisin E variants 

(88). The mutations at base positions 1 1 07 in RC 1 and 995 in RC2 (Figure 14), giving 
rise to amino acid substitutions Asn218/Ser (N218S) and Asnl81/Asp (N181 ID), lead to 
improvements in subtiiisin E thermostability; the remaining mutations, both synonymous 
and nonsvnonymous. have no detectable effects on thermostability. At 65°C, the single 
10 variants NISID and N218S have approximately 3-fold and 2-fold longer half-lives, 

respectively, than wild subtiiisin E. and variants containing both mutations have half-lives 
that are 8-fold longer (88V The different half-lives in a population of subtiiisin E variants 
can therefore be used to estimate the efficiency by which sequence information is 
combined. In particular, recombination between these two mutations (in the absence of 
15 point mutations affecting thermostability) should generate a library in which 25% of the 
population exhibits the thermos/ability of the double mutant. Similarly, 25% of the 
population should exhibit wild-type like stability, as N181D andN218S are eliminated at 
equal frequency. We used the fractions of the recombined population as a diagnostic 

20 A. Methods 

The strategy underlying this example is shown in Fig. 15. 
Subtiiisin E thermostable mutant genes RC1 and RC2 (Fig. 14) are 986-bp 
fragments including 45 nt of subtiiisin E prosequence. the entire mature sequence and 113 
nt after the stop codon. The genes were cloned between Bam HI and Nde I in E. colilB. 
25 subtilis shunle vector pBE3. resulting in pBE3-l and pBE3-2. respectively. Plasmid 
DNA pBE3-l and pBE3-2 was isolated from E.coli SCSI 10. 

About 5.0 ug of ummethyiated pBE3-l and pBE3-2 DNA were digested 
with Bam HI and Sdc I. respectively, at 37°C for 1 hour. After agarose gel separation, 
equimoiar concentrations (2.0nM) of the linearized unmethylated pBE3-l and pBE3-2 
30 were mixed in 1 x SSPE buffer ( 1 80 mM NaCl. 1 .0 mM EDTA. 1 0 mM NaH.PCu. pH 

7.4 i. After heating at 96 = C for 10 minutes, the reaction mixture was immediately cooled 
at 0°C for 5 min. The mixture was incubated at 68°C for 2 hr for heterodupiexes to form. 
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One microliter of the reaction mixture was used to transform 50 jil of £ 
coli ES 1301 mutS. E. coli SCSI 10 and £ co//HB101 competent cells. 

The transformation efficiency with £ coli HB101 competent cells was 
about ten times higher than that of E.coli SCSI 10 and 15 times higher than that of £ coli 
5 ESI 301 mutS. But in all these cases, the transformation efficiencies were 10-250 times 
lower than that of the transformation with closed, covalent and circular control pUC19 
plasmids. 

Five clones from £ coli SCSI 10 mutant library and five from £ coli 
ESI 301 mutS library were randomly chosen, and plasmid DNA was isolated using a 
1 0 QIAprep spin plasmid miniprep kit for further DNA sequencing analysis. 

About 2.000 random clones from £ coli HB101 mutant library were 
pooled and total plasmid DNA was isolated using a Q1AGEN-100 column. 0.5-4.0 ug of 
the isolated plasmid was used to transform Bacillus subtil is DB428 as described 
previously (88). 

1 5 About 400 transformants from the Bacillus subtilis DB428 library were 

subjected to screening. Screening was performed using the assay described previously 
(88), on succinyl-AIa-Ala-Pro-Phe-p-nitroanilidc. B. subtilis DB428 containing the 
plasmid library were grown on LB plates containing kanamycin (20 jag/ml) plates. After 
1 8 hours at 37°C single colonies were picked into 96-weil plates containing 200 ill 
20 SG/kanamycin medium per well. These plates were incubated with shaking at 37°C for 
24 hours to let the cells to grow to saturation. The cells were spun down, and the 
supernatants were sampled for the thermostability assay. 

Two replicates of 96-well assay plates were prepared for each growth plate 
by transferring 10 ul of supernatant into the replica plates. The subtilisin activities were 
25 then measured by adding 100 [i\ of activity assay solution (0.2 mM succinyl-Ala-Ala-Pro- 
Phe-p-mtroanilide. 100 mM Tris-HCl, 10 mM CaCK pH 8.0. 37°C). Reaction velocities 
were measured at 405 nm to over 1.0 min in a ThermoMax microplate reader (Molecular 
Devices. Sunnyvale CA). Activity measured at room temperature was used to calculate 
the fraction of active clones (clones with activity less than 10% of that of wild type were 
30 scored as inactive). Initial activity (A,) was measured after incubating one assay piate at 
65°C for 10 minutes by immediately adding 100 ul of prewarmed (37°C) assay solution 
(0.2 mM succmyl-Aia-Ala-Pro-Phe-p-nitroamlide. I 00 mM Tris-HCl. pH 8.0. 10 mM 
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CaCl : , pH 8.0) into each well. Residual activity {An was measured after 40 minute 
incubation. 

B. Results 

In vuro heterodupiex formation and in vivo repair was earned out as 
5 described above. Five clones from E, coii SCSI 10 mutant library and five from E. coli 
ESI 301 mutS libraries were selected at random and sequenced. Fig. 14 shows that four 
out of the ten clones were different from the parent genes. The frequency of occurrence 
of a particular point mutation from parent RC1 or RC2 in the resulting genes ranged from 
0% to 50%, and the ten point mutations in the heterodupiex have been repaired without 
1 0 strong strand-specific preference. 

Since none of the ten mutations locates within the dem sue. the mismatch 
repair appears generally done via the E. coii long-patch mismatch repair systems. The 
system repairs different mismatches in a strand-specific manner using the state of N6- 
methvlation of adenine in GATC sequences as the major mechanism for determining the 
1 5 strand to be repaired. With heteroduplexes methylated at GATC sequences on only one 
DNA strand, repair was shown to be highly biased to the unmethylated strand, with the 
methylated strand serving as the template for correction. If neither strand was 
methylated, mismatch repair occurred, but showed little strand preference (23, 24). These 
results shows that it is preferable to demethylate the DNA to be recombined to promote 
20 efficient and random repair of the heteroduplexes. 

The rates of subtihsin E thermo-inactivation at 65°C were estimated by 
analyzing the 400 random clones from the Bacillus subtilis DB42S library. The 
thermostabilities obtained from one 96-weli plate are shown in Figure 16. plotted in 
descending order. About 12.9% of the clones exhibited thermostability comparable to the 
25 mutant with the NI81D and N2 1 8S double mutations. Since this rate is only half of that 
expected for random recombination of these two markers, it indicates that the two 
mismatches at positions 995 and 1 107 within the heteroduplexes have been repaired with 
lower position randomness. 

Sequence analysis of the clone exhibiting the highest thermostability 
30 amone the screened 400 transformants trom the E. coli SCSI 10 heterodupiex library 
confirmed the presence of both N181D and K21SS mutations. Among the 400 
transformants from the B. subtilis DB428 library that were screened, approximately 91% 
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of the clones expressed N181D- and/or N218S-type enzyme stabilities, while about 8.0% 
of the transformants showed only wild-type subtihsin E stability. 

Less than 1 .0% inactive clone was found, indicating that few new point 
mutations were introduced in the recombination process. This is consistent with the fact 
5 that no new point mutations were identified in the ten sequenced genes (Figure 14). 
While point mutations may provide useful diversity for some in vitro evolution 
applications, they can also be problematic for recombination of beneficial mutations, 
especially when the mutation rate is high. 

10 EXAMPLE 4. Optimizing Conditions For The Heteroduplex Recombination. 

We have found that the efficiency of heteroduplex recombination can 
differ considerably from gene to gene [17,57]. In this example, we investigate and 
optimize a variety of parameters that improve recombination efficiency. 
DNA substrates used in this example were site-directed mutants of green fluorescent 

1 5 protein from Aequorea victoria. The GFP mutants had a stop codon(s) introduced at 

different locations along the sequence that abolished their fluorescence. Fluorescent wild 
type protein could be only restored by recombination between two or more mutations. 
Fraction of fluorescent colonies was used as a measure of recombination efficiency. 

20 A. Methods 

About 2-4 jj.g of each parent plasmid was used in one recombination 
experiment. One parent plasmid was digested with Pst I endonuclease another parent 
with EcoRl. Linearized plasmids were mixed together and 20 x SSPE buffer was added to 
the final concentration lx (180 mM NaCL 1 mM EDTA. 10 mM NaH 2 P0 4 , pH 7.4). The 
25 reaction mixture was heated at 96°C for 4 minutes, immediately transferred on ice for 4 
minutes and the incubation was continued for 2 hours at 68°C. 

Target genes were amplified in a PCR reaction with primers corresponding 
to the vector sequence of pGFP plasmid. Forward primer: 5*- 
CCGACTGGAAAGCGGGCAGTG-3\ reverse pnmer 5*- 
30 CGGGGCTGGCTTAACTATGCGG-3 * . PCR products were mixed together and purified 
using Qiagen PCR purification kit. Purified products were mixed with 20 x SSPE buffer 
and hybridized as described above. Annealed products were precipitated with ethanol or 
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purified on Qiagen columns and digested with EcoRl and Psil enzymes. Digested 
products were ligated into Pst\ and EcoRl digested pGFP vector. 

dUTP was added into PCR reaction at final concentrations 200 uM. 40uM 
8 u\l. 1.6 uM. 0.32 uM. PCR reaction and subsequent cloning procedures were 
performed as described above. 

Recombinant plasmids were transformed into XL10 E. coii strain by a 
modified chemical transformation method. Cells were plated on ampicillin containing LB 
acar plates and grown overnight at 37°C, followed by incubation at room temperature or 
at 4°C until fluorescence developed. 

B. Results . 

1 . Effect of ligation on recombination efficiency . 

Two experiments have been performed to test the effect of breaks in the 
DNA heteroduplex on the efficiency of recombination. In one experiment heteroduplex 
plasmid was treated with DNA ligase to close all existing single-strand breaks and was 
transformed in identical conditions as an unligated sample (see Table 1). The ligated 
samples show up to 7-fold improvement in recombination efficiency over unligated 
samples. 

In another experiment, dUTP was added into PCR reaction to introduce 
0 additional breaks into DNA upon repair by uracyl N-glycosylase in the host cells. Table 
2 shows that dUMP incorporation significantly suppressed recombination, the extent of 
suppression increasing with increased dUTP concentration. 

2. Effect of plasmid size on the efficiency of 
5 heteroduplex formation . 

Plasmid size was a significant factor affecting recombination efficiency. 
Two plasmids pGFP (3.3 kb) and a Bacillus shuttle vector pCTl (about 9 kb> were used 
in preparing circular heteroduplex-like plasmids following traditional heteroduplex 
protocol. For the purpose of this experiment (to study the effect of plasmid size on 
mj duplex formation), both parents had the same sequences. While pGFP formed about 30 
40% of circular plasmid. the shuttle vector yielded less than 10% of this form. 
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Increase in plasmid size decreases concentration of the ends in the vicinity 
of each and makes annealing of very long (>0.8 kb) ends that are single-stranded more 
difficult. This difficulty is avoided by the procedure shown in Fig. 3, in which 
heteroduplex formation occurs between substrates in vector-free form, and, 
heteroduplexes are subsequently inserted into a vector. 



3. Efficiency of Recombination vs. Distance Between Mutations 
A series of GFP variants was recombined pairvvise to study the effect of 
distance between mutations on the efficiency of recombination. Parental genes were 
10 amplified by PCR. annealed and ligated back into pGFP vector. Heteroduplexes were 
transformed into XL 10 E.coii strain. 

The first three columns in Table 3 show the results of three independent 
experiments and demonstrate the dependence of recombination efficiency on the distance 
between mutations. As expected recombination becomes less and less efficient for very 
1 5 close mutations. 

However, it is still remarkable that long-patch repair has been able to 
recombine mutations separated by only 27 bp. 

The last line in Table 3 represents recombination between one single and 
one double mutants. Wild type GFP could only be restored in the event of double 
20 crossover with each individual crossover occurring in the distance of 99 bp only, 

demonstrating the ability of this method to recombine multiple, closely-spaced mutations. 

4. Elimination Of The Parental Double Strands 
From Heteroduplex Preparations. 

25 Annealing of substrates in vector-free form offers size-advantages relative 

to annealing of substrates as components of vectors, but does not allow selection for 
heteroduplexes relative to homoduplexes simply by transformation into host. 
Asymmetric PCR reactions with only one primer for each parent seeded with appropriate 
amount of previously amplified and purified gene fragment were run for 100 cycles. 

30 ensuring a 100-fold excess of one strand over another. Products of these asymmetrical 
reactions were mixed and annealed together producing only a minor amount ot 
nonrecombinant duplexes. The last column in Table 3 shows the recombination 
efficiency obtained from these enriched heteroduplexes. Comparison of the first three 
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columns with the fourth one demonstrates the improvement achieved by asymmetric 
synthesis of the parental strands. 

While the foregoing invention has been described in some detail for 
purposes of clarity and understanding, it will be ciear to one skilled in the an from a 
reading of this disclosure that various changes in form and detail can be made without 
departing from the true scope of the invention. All publications and patent documents 
cited in this application are incorporated by reference in their entirety for all purposes to 
the same extent as if each individual publication or patent document were so individually 
denoted. 
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WHAT IS CLAIMED IS: 

I. A method for evolving a polynucleotide toward acquisition of a 
desired property, comprising 



3 (a) incubating a population of parental polynucleotide variants under 

4 conditions to generate annealed polynucleotides comprising heteroduplexes; 

5 (b) exposing the heteroduplexes to a cellular DNA repair system to 

6 convert the heteroduplexes to parental polynucleotide variants or recombined 

7 polynucleotide variants: 

8 (c) screening or selecting the recombined polynucleotide variants for the 

9 desired property. 

1 2. The method of claim i, wherein the heteroduplexes are exposed to 

2 the cellular DNA repair system in vitro. 

1 3. The method of claim 2. wherein the cellular DNA repair system 

2 comprises cellular extracts. 

1 4. The method of claim 1, further comprising introducing the 

2 heteroduplexes into cells, whereby the heteroduplexes are exposed to the DNA repair 

3 system of the cells in vivo. 

1 5. The method of claim 4. wherein the annealed polynucleotides 

2 further comprise homoduplexes and the introducing step selects for transformed cells 

3 comprising the heteroduplexes relative to transformed cells comprising homoduplexes. 

1 6. The method of claim 4, wherein a first polynucleotide variant is 

2 provided as a component of a first vector, and a second polynucleotide variant is provided 
as a component of a second vector, and the method further comprises convening the first 

4 and second vectors to linearized forms in which the first and second polynucleotide 

5 variants occur at opposite ends, whereby in the incubating step single-stranded forms of 
the firs: iineanzed vector reanneai with each other to form iinear first vector, single- 
stranded forms of the second linearized vector reanneai with each other to form iinear 

S second vector, and singie-stranded Iineanzed forms of the first and second vectors anneal 

c > with each to form a circular heteroduplex bearing a nick in each strand, and the 
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1 introducing step selects for transformed cells comprising the circular 

2 heteroduplexes relative to the linear first and second vector. 

1 7. The method of claim 6, wherein the first and second vectors are 

2 converted to linearized forms by PCR. 

1 8. The method of claim 6, wherein the first and second vectors are 

2 converted to linearized forms by digestion with first and second restriction enzymes. 

1 9. The method of claim L wherein the population of polynucleotide 

2 variants are provided in double stranded form, and the method further comprising 

3 converting the double stranded polynucleotides to single stranded polynucleotides before 

4 the annealing step. 

1 10. The method of claim 1, wherein the converting step comprises: 

2 conducting asymmetric amplification of the first and second double 

3 stranded polynucleotide variants to amplify a first strand of the first polynucleotide 

4 variant, and a second strand of the second polynucleotide variant, whereby the first and 

5 second strands anneal in the incubating step to form a heteroduplex. 

1 11. The method of claim 10. wherein the first and second double- 

2 stranded polyncueltoides variants are provided in vector-free form, and the method 

3 further comprises incorporating the heteroduplex into a vector. 

1 12. The method of claim 4 wherein the population of polynucleotides 

2 comprises first and second polynucleotides provided in double stranded form, and the 

3 method further comprises incorporating the first and second polynucleotides as 

4 components of first and second vectors, whereby the first and second polynucleotides 

5 occupy opposite ends of the first and second vectors, whereby in the incubating step 

6 single-stranded forms of the first linearized vector reanneal with each other to form linear 

7 first vector, single-stranded forms of the second linearized vector reanneal with each 

8 other to form linear second vector, and single-stranded linearized forms of the first and 

9 second vectors anneal with each to form a circular heteroduplex bearing a nick in each 
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1 strand, and the introducing step selects for transformed celis comprises the 

2 circular heterodupiexes relative to the linear first and second vector. 

1 13. The method of claim 4. further comprising sealing nicks in the 

2 heterodupiexes to form covalently-closed circular heterodupiexes before the introducing 

3 step. 

1 14. The method of claim 1 1 , wherein the first and second 

2 polynucleotides are obained from chromosomal DNA.. 



1 15. The method of claim 1. further comprising repeating steps (a)-(c) 

2 whereby the incubating step in a subsequent cycle is performed on recombinant variants 

3 from a previous cycle. 

1 16. The method of claim 1. wherein the polynucleotide variants encode 

2 a polypeptide. 

1 17. The method of claim 1 . wherein the population of polynucleotide 

2 variants comprises at least 20 variants. 

1 18. The method of claim 1 , wherein the population of polynucleotide 

2 variants are at least 10 kb in length. 

1 19. The method of claim 1. wherein the population of polynucleotide 

2 variants comprises natural variants. 

1 20. The method of claim 1, wherein the population of polynucleotides 

2 comprises variants generated by mutagenic PCR. 

1 21. The method of claim 1 . wherein the population of polynucleotide 

2 variants comnses variants generated by site directed mutagenesis. 

] 22. The method of claim 1. wherein the ceils are bactenai celis. 

1 23. The method of claim i. further comprising at least partially 

2 demethylating the population of variant polynucleotides. 
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1 24. The method of claim 23. whether the at least partially 

2 demethylating step is performed by PCR amplification of the population of variant 

3 polynucleotides. 

1 25. The method of claim 23. wherein the at least partially 

2 demethylating step is performed by amplification of the population of variant 

3 polynucelotides in host cells. 

1 26. The method of claim 25, wherein the host cells are defective in a 

2 gene encoding a methylase enzyme. 

1 27. The method of claim L wherein the population of variant 

2 polynucleotide variants comprises at least 5 polynucleotides having at least 90% sequence 

3 identity with one another. 

1 28. The method of claim 1, further comprising isolating a screened 

2 recombinant variant. 

1 29. The method of claim 28, further comprising expressing a screened 

2 recombinant variant to produce a recombinant protein. 

1 30. The method of claim 29 further comprising formulating the 

2 recombinant protein with a earner to form a pharmaceutical composition. 

1 31. The method of claim 1 , wherein the polynucleotide variants encode 

2 enzymes selected from the group consisting of proteases, lipases, amylases, cutinases, 

3 cellulases, amylases, oxidases, peroxidases and phytases. 

1 32. The method of claim 1. wherein the polynucleotide variants encode 

2 a polypeptide selected from the group consisting of insulin. ACTH. glucagon. 

3 somatostatin, somatotropin, thymosin, parathyroid hormone, pigmentary hormones. 

4 somatomedin, enhropoietin. luteinizing hormone, chorionic gonadotropin, hyperthalmic 

5 releasing factors, antidiuretic hormones, thyroid stimulating hormone, reiaxm. interferon. 

6 thrornbopoietin (TPCM. and prolactin. 



SNSDOCID -:WO 



OQ29902A1 I > 



WO 99/29902 



PCT /US98/25698 



-47- 

33. The method of claim I. wherein the polynucleotide variants 
encode a plurality of enzymes forming a metabolic pathway. 

34. The method of claim 1. wherein the polynucleotide variants are in 
concatemenc form. 
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C7GCAGCCTGCCGACC7G7TCG7GGTGGTGATCGCGGCCGCGCTGGCCGCC37CGCGGTC 
GCCGCCGCCGGGCCGATCGAGTTGGTCGCCTTCGTCSTGCCGCAGATCGCC^TGCGGCTC 
TGCGGCGGCAGCCGGCCGCCCCTGCTCGCCTCGGCGATGCTCGGCGCGCTGCTCGTGGTC 
GGCGCCGACCTGGTCGCTCAGATCGTGGTGGCGCCGAAGGAGCTGCCGGTCGGCCTGCTC 
ACCGCGATGATCGGCAGCCCGTACCTGCTCTGGCTCCTGCTTCGGCGA7CAAGAAAGGTG 
AGCGGATGAACGCCCGCCTGCG7GGCGAGGGCCTGCACCTCGCGTACGGGGACCTGACCG 

3 £ 1 TGATCGACGGCCTCGACGTCG ACGTGCACGACGGGCTGGTCACCACCATCATCGGGCCCA 
421 ACGGGTGCGGCAAGTCGACGCTGCTCAAGGCGCTCGGCCGGCTGCTGCGCGCGACCGGCG 

4 81 GGCAGGTGCTGCTGGACGGCCGCCGCATCGACCGGACCCCCACCCGTGAGGTGGCCCGGG 

5 41 TGCTCGGCGTGCTGCCGCAGTCGCCCACCGCGCCCGAAGGGCTCACCGTCGCCGACCTGG 

6 C 1 TGA7GCGCGGCGGGGACCCGCACCAGACC7GGTTCCGGCAGTGGTCGCGGGACGACGAGG 

6 £ 1 ACCAGGTCGCCGACGCGGTGCGCTGGACCGACATGCTGGCGTACGCGGACCGCCCGGTGG 

7 21 ACGCCCTCTCCGGCGGTC AGCGCCAGCGCGCCTGGATCAGC ATGGCGCTGGCCCAGGGCA 
7 6 1 CCGACCTGCTGCTGGTGGACGAGCCGACCACCTTCCTCGACCTGGCCCACCAGATCGACG 
B 4 1 TGCTGGACCTGGXCCCCCGGCTGCACGCCGAGATGGGCCGGACCGTGGTGATGGTGCTGC 
9 01 ACGACCTGAGCCTGGCCGCCCGGTACGCCGACCGGCTGATCGCGATGAAGGACGGCCGGA 

9 61 TCGTGGCGAGCGGGGCGGGGGACGAGG7GCTCACCCCGGCGCTGCTGGAGTCGGTCTTCG 
1C21 GGCTGCGCGCGATGGTGGTGCCCGACCCGGCGACCGGCACCCCGCTGGTGATCCCCCTGC 

10 81 CGCGCCCCGCCAGCTCGG7GCGGGGCTGAAATCGATGAGCGTGGTTGCTTC ATCGGCCTG 
1141 CCGAGCGATGAGAGTATGTGGGCGGTAGAGCGAGTCTCGAGGGGGAGATGCCGCCGTGAC 

V T 

1201 G7CCTCGTACATGCGCC7GAAAGCAGCAGCGATCGCC7TCGG7GTGATCG7GGCGACCGC 
5 SSYMRLKAAAIAFGVIVATA 

12 61 AGCCGTGCCGTCACCCGCTTCCGGCAGGGAACATGACGGCG3CTA7GCGGG CCTGATCCG 
23 AVPSPASGREHDGGYAALIR 

13 21 CCGGGCC7CGTACGGCGTCCCGGACATCACCGCCGACGACT7CGGGAGCC7CGGTTTCGG 
4 3 R A S Y G V ? K Z 7 A D 2 F G S L G F 3 



.441 CAAGGGTGAGCGG7CGG: 



;TGCGACGGG-GCCGGACGACGCCGA7G7GGGCAG 
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83 NGERSRWFGATGPDDADVRS 

15 01 CGACCTCTTCCACCGCAAGGCGATCGACGACCGCGTCGCCGAGCGGCTCCTCGAAGGGCC 
103 DLFHRKAIDDRVAERLLEGF 

15 61 CCGCGACGGCGTGCGGGCGCCGTCGGACGACGTCCGGGACCAGATGCGCGGCTTCGTCGC 
123 RDGVRAP SDDVRDQMRG FVA 

1621 CGGCTACAACCACTTCCTACGCCGCACCGGCGTGCACCGCCTGACCGACCCGGCGTGCCG 
143 GYNHFLRRTGVHRLTDPACR 

i 6 B 1 CGGCAAGGCCTGGGTGCGCCCGCTCTCCGAGATCGATCTCTGGCGTACGTCGTGGGACAG 
163 GKAWVRPLSEIDLWRTSWDS 

17 41 CATGGTCCGGGCCGGTTCCGGGGCGCTGCTCGACGGCATCGTCGCCGCGACGCCACCTAC 
183 MVRAG5GALLDGTVAATP P T 

16 01 AGCCGCCGGGCCCGCGTCAGCCCCGGAGGCACCCGACGCCGCCGCGATCGCCGCCGCCCT 
203 AAGPASAP EAPDAAAIAAAL 

18 61 CGACGGGACGAGCGCGGGC ATCGGCAGCAACGCGTACGGCCTCGGCGCGCAGGCCACCGT 
223 DGTSAG1GSNAYGLGAQATV 

IS 21 GAACGGCAGCGGGATGGTGCTGGCCAACCCGCACTTCCCGTGGCAGGGCGCCGCACGCTT 
243 NGSGKVLANPHFPWQGAARF 

19 El CTACCGGATGCACCTCAAGGTGCCCGGCCGCTACGACGTCGAGGGCGCGGCGCTGATCGG 
2 63 YRKHLKVFGRYDVEGAALIG 

2 0 41 CGACCCGATCATCGGGATCGGGCACAACCGC ACGGTCGCCTGGAGCCACACC3TCTCCAC 

2 S3 D P I I G I G H N R 7 V A W S K T V S 7 

210 1 CGCCCGCCGGTTCGTGTGGC ACCGCC7GAGCC7GGTGCCCGGCGACCCCACC7CCTATTA 

3 03 ARRFVWHRLSLVPGDPTSYY 

Fig. 13B 
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216 1 CC~CCACGGCCGGCCCGAGCGGATGCGCGCCCGCACGGTCACGGTCCAGACCGGCA3CGC 
32 3 V D G R P HI R M R A R T V T V Q T G S G 




22 El CTTCGACTGGACGCCGGCCACCG2GTACGGCATCACCGACGTCAACGCGGGCAACAACCG 
363 7 D W T P A T A Y A I T D V N A G N N R 

2 3^1 CGCCTTCGACGGGTGGCTGCGGATGGGCCAGGCCAACGACGTCCGGGCGCTCAAGGCGGT 
j£J A F D G W L R K G Q A D V R A L K A V 

2 4 01 CCTCGACCGGCACCAGTTCCTGCCCTGGGTC AACGTGATCGCCGCCGACGCGCGCCGCGA 
LDRHQFI-PWVNVIAADARGE 

2 4 61 GGCCCTCTACGGCGATCATTCGGTCGTCCCCCGGGTGACCGGCGCGCTCGCTGCCGCCTG 
4Z3 ALYGDHSVVFRVTGALAAAC 




2 5B1 CCGGTCGGACTGCGCGCTCGGCGCCGACCCCGACGCCGCGGTCCCGGGGATTCTCGGCCC: 
4 5 j HSDCALGADPDAAVPGI^GP 

2 6 41 GGCGAGCCTGGCGGTGCGGTTCCGCGACGACTACGTCACCAACTCCAACGACA3TCAC7G 

4 53 ASLPVRPRDDYVTNSNDSHW 

2 7 01 GCTGGCCAGCCCGGCCGGGCCGCTGGAAGGCTTCCCGCGGATCCTCGGCAACGAACGCAC 

5 03 i.ASPAAPLEGFPRILGNERT 



27 61 ccggcgcagcctgcg: 



VTCCAGCAGCGCCTCGCCGGCAC 



597 p R S l F . r R L G L C 0 I C C R 



A G 
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2 8 B 1 CCGGATGCACGGCGCCGAACTCGCCCGCGACGACCTGGTCGCGCTCTGCCGCCGCCAGCC 
563 RMHGAELARDDLVALCKRQP 

29 41 G ACCGCGACCGCCTCGAACGGCGCGATCGTCGACCTCACCGCGGCCTCCACGGCGCTGTC 
593 TATASNGAIVDLTAACTALS 

3 0 0 1 CCGCTTCGATGAGCGTGCCGACCTGGACAGCCGGGGCGCGCACCTGTTCACCGAGTTCGC 
603 R F D E RAD L D S R G _ A H LFT E FA 

30 61 CCTCGCGGGCGGAATCAGGTTCGCCGACACCTTCGAiGGTGACCGATCCGGTACGCACCCC 
623 LAGGI RFADTFE VTDPVRTP 

3121 GCGCCGTCTGAACACCACGGATCCGCGGGTACGGACGGCGCTCGCCGACGCCGTGCAACG 
643 RRLNTTDP RVRTALADAVQR 

3181 GCTCGCCGGCATCCCCCTCGACGCGAAGCTGGGAGACATCCACACCGACAGCCGCGGCGA 
663 LAGI PL-DAKLGDI HTDS RGE 

324 1 ACGGCGCATCCCCATCCACGGTGGCCGCGGGGAAGCAGGCACCTTCAACGTGATCACCAA 
683 RRIPIHGGRGEAGTFNVITN 

3 3 01 CCCGCTCGTGCCGGGCGTGGGATACCCGCAGGTCGTCCACGGAAC ATCGTTCGTGATGGC 
703 PLVPGVGYPQVVHGTSFVMA 

3 3 61 CGTCGAACTCGGCCCGCACGGCCCGTCGGGACGGCAGATCCTCACCTATGCGCAGTCGAC 
723 VELGPHGPSGRQI L T Y A Q S T 

3 4 21 GAACCCGAACTCACCCTGGTACGCCGACCAGACCGTGCTCTACTCGCGGAAGGGCTGGGA 
742 NPNSPWYADQTVLVSRKGWD 

3 4 81 CACCATCAAGT ACACCGAG^CGCAGATCGCGGCCGACCCGAACCTGCGCGTCTACCGGGT 
7 63 T I K Y T E A Q r A A D P N L R V Y R V 

3 54 1 GGCACAGCGGGGACGCTGACCCACGTCACGCCGGCTCGGCCCGTGCGGGGGCGCAGGGCG 

Fig. 13D 
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7B3 A Q R G K 

3 € C 1 CCGATCGTCTCTGCATCGCCGGTCAGCCGGGGCCTGCGTCGACCGGCGGCGGCCGG7CGA 

3 6 € I CGCCCGCGTCCCGGCGCAGCGACTGGCTGAAGCGCCAGGCGTCGGCGGCCCGGGGCAGGT 

37 21 7G77GAACA7CAGG7ACGCCGGGCCGCCG7CGAGGA7GCCGGCGAGGTGTGCCAGC7CGG 

3 7 81 C A7CCG7G7ACACATGCCGGGCGCCGGTGATGCCG7GCAGCCGG7AATAGGCCATCGGCG 

3 6 41 TCAGACTGCGGCGCAGGAACGGGTCGGCGGCGTGGGTCAGGTCCAGCTCCTGGCACAAGC 

3 901 CGTCGACCACCTCG7CCGGCCACGGGCCGCGCGGCTCCCACAACAGCCGGACACCGGCCG 

3 9 61 GCCGGCGCGGTCGGGCGCAG AACTCACGCAGTCGCGCGATGGCGGG7TCGGTCGGCCGGA 

4 0 21 AAC7CGCCGGGCAC7GCAG 



Fig. 13E 
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METHOD FOR CREATING POLYNUCLEOTIDE AND 
POLYPEPTIDE SEQUENCES 

CROSS-REFERENCES TO RELATED APPLICATIONS 
This application derives priority from USSN 60/067908, filed 
December 8, 1997, which is incorporated by reference in its entirety for all purposes. 

) TECHNICAL FIELD 

The invention resides in the technical field of genetics, and more 
specifically, forced molecular evolution of polynucleotides to acquire desired properties. 

BACKGROUND 

5 A variety of approaches, including rational design and directed evolution, 

have been used to optimize protein functions (1,2). The choice of approach for a given 
optimization problem depends, in part, on the degree of understanding of the relationships 
between sequence, structure and function. Rational redesign typically requires extensive 
knowledge of a structure-function relationship. Directed evolution requires little or no 
0 specific knowledge about structure-function relationship; rather, the essential features is a 
means to evaluate the function to be optimized. Directed evolution involves the 
generation of libraries of mutant molecules followed by selection or screening for the 
desired function. Gene products which show improvement with respect to the desired 
property or set of properties are identified by selection or screening. The gene(s) 
25 encoding those products can be subjected to further cycles of the process in order to 

accumulate beneficial mutations. This evolution can involve few or many generations, 
depending on how far one wishes to progress and the effects of mutations typically 
observed in each generation. Such approaches have been used to create novel functional 
nucleic acids (3, 4), peptides and other small molecules (3 k antibodies (3 >. as well as 
30 enzymes and other proteins (5, 6. " ). These procedures are fairly tolerant to inaccuracies 
and noise in the function evaluation ('7}. 
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Several publications have discussed the role of gene recombination in 
directed evolution (see WO 97/07205. WO 98/42727, US 5807723. US 5.721.367, 
US 5.776.744 and WO 98/41645 US 5.81 1.238. WO 98/41622. WO 98/41623. and 
US 5.093.257). 

A PCR-based group of recombination methods consists of DNA shuffling 
[5, 6], staggered extension process [89. 90] and random-priming recombination [87]. 
Such methods typically involve synthesis of significant amounts of DNA during 
assembly/recombination step and subsequent amplification of the final products and the 
efficiency of amplification decreases with gene size increase. 

Yeast cells, which possess an active system for homologous 
recombination, have been used for in vivo recombination. Cells transformed with a 
vector and partially overlapping inserts efficiently join the inserts together in the regions 
of homology and restore a functional, covalently-closed plasmid [91]. This method does 
not require PCR amplification at any stage of recombination and therefore is free from the 
size considerations inherent in this method. However, the number of crossovers 
introduced in one recombination event is limited by the efficiency of transformation of 
one cell with multiple inserts. Other in vivo recombination methods entail recombination 
between two parental genes cloned on the same plasmid in a tandem orientation. One 
method relies on homologous recombination machinery of bacterial cells to produce 
chimeric genes [92]. A first gene in the tandem provides the N-terminal part of the target 
protein, and a second provides the C-terminal part. However, only one crossover can be 
Generated by this approach. Another in vivo recombination method uses the same tandem 
organization of substrates in a vector [93]. Before transformation into E. coli cells, 
plasmids are linearized by endonuclease digestion between the parental sequences. 
5 Recombination is performed in vivo by the enzymes responsible for double-strand break 
repair. The ends of linear molecules are degraded by a 5'->3' exonuclease activity, 
followed by annealing of complementary single-strand 3* ends and restoration of the 
double-strand plasmid [94]. This method has similar advantages and disadvantages of 
tandem recombination on circular plasmid. 

0 

SUMMARY OF THE INVENTION 
The invention provides methods for evolving a polynucleotide toward 
acquisition of a desired property. Such methods entail incubating a population of 
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parental polynucleotide variants under conditions to generate annealed polynucleotides 
comprises heterodupiexes. The heterodupiexes are then exposed to a cellular DNA 
repair system to conven the heterodupiexes to parental polynucleotide variants or 
recombined polynucleotide variants. The resulting polynucleotides are then screened or 

5 selected for the desired property. 

In some methods, the heterodupiexes are exposed to a DNA repair system 
in vitro. A suitable repair system can be prepared in the form of cellular extracts. 

In other methods, the products of annealing including heterodupiexes are 
introduced into host cells. The heterodupiexes are thus exposed to the host cells' DNA 

0 repair system in vivo. 

In several methods, the introduction of annealed products into host cells 
selects for heterodupiexes relative to transformed cells comprising homoduplexes. Such 
can be achieved, for example, by providing a first polynucleotide variant as a component 
of a first vector, and a second polynucleotide variant is provided as a component of a 
5 second vector. The first and second vectors are converted to linearized forms in which 
the first and second polynucleotide variants occur at opposite ends. In the incubating 
step, single-stranded forms of the first linearized vector reanneal with each other to form 
linear first vector, single-stranded forms of the second linearized vector reanneal with 
each other to form linear second vector, and single-stranded linearized forms of the first 
20 and second vectors anneal with each to form a circular heteroduplex bearing a nick in 

each strand. Introduction of the products into cells thus selects for cirular heterodupiexes 
relative to the linear first and second vector. Optionally, in the above methods, the first 
and second vectors can be converted to linearized forms by PCR. Alternatively, the first 
and second vectors can be converted to linearized forms by digestion with first and 

second restriction enzymes. 

In some methods, polynucleotide variants are provided in double stranded 
form and are converted to single stranded form before the annealing step. Optionally, 
such conversion is by conducting asymmetric amplification of the first and second 
double stranded polynucleotide variants to amplify a first strand of the first 
$0 polynucleotide variant, and a second strand of the second polynucleotide variant. The 
first and second strands anneal in the incubating step to form a heteroduplex. 

in some methods, a population of polynucleotides comprising firs; and 
second polynucleotides is provided in double stranded form, and the method further 
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compnses incorporating the first and second polynucleotides as components of first and 
second vectors, whereby the first and second polynucleotides occupy opposite ends of 
the first and second vectors. In the incubating step single-stranded forms of the first 
linearized vector reanneal with each other to form linear first vector, single-stranded 
forms of the second linearized vector reanneal with each other to form linear second 
vector, and single-stranded linearized forms of the first and second vectors anneal with 
each to form a circular heteroduplex bearing a nick in each strand. In the introducing 
step selects for transformed cells comprises the circular heteroduplexes relative to the 
linear first and second vector. 

In some methods, the first and second polynucleotides are obtained from 
chromosomal DNA. In some methods, the polynucleotide variants encode variants of a 
polypeptide. In some methods, the population of polynucleotide variants comprises at 
least 20 variants. In some methods, the population of polynucleotide variants are at least 
1 0 kb in length. 

In some methods, the polynucleotide variants comprises natural variants. 
In other methods, the polynucleotide variants comprise variants generated by mutagenic 
PCR or cassette mutagenesis. In some methods, the host cells into which heteroduplexes 
are introduced are bactciial cells. In some methods, the population of variant 
polynucleotide variants comprises at least 5 polynucleotides having at least 90% sequence 
identity with one another. 

Some methods further comprise a step of at least partially demethylating 
variant polynucleotides. Demethylation can be achieved by PCR amplification or by 
passaging variants through methylation-deficient host cells. 

Some methods include a further step of sealing one or more nicks in 
5 heteroduplex molecules before exposing the heteroduplexes to a DNA repair system. 
Nicks can be sealed by treatment with DNA hgase. 

Some methods further comprise a step of isolating a screened recombinant 
polynucleotide ariant. In some methods, the polynucleotide variant is screened to 
produce a recombinant protein or a secondary metabolite whose production is catalyzed 
0 thereby. 

In some methods, the recombinant protein or secondary metabolite is 
formulated with a earner to form a pharmaceutical composition. 
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Figure 9 shows the physical maps of Actinoplanes uiahensis ECB 
deacylase mutants with enhanced specific activity ((a) is pM7-2 for Mutant 7-2. and (b) is 
pM16 for Mutant 16). 

Figure 10 illustrates the process used for Example 2 to recombine 
5 mutations in Mutant 7-2 and Mutant 16 to yield ECB deacylase recombinant with more 
enhanced specific activity. 

Figure 1 1 shows specific activities of wild-type ECB deacylase and 
improved mutants Mutant 7-2, Mutant 16 and recombined Mutant 15. 

Figure 12 shows positions of DNA base changes and amino acid 
1 0 substitutions in recombined ECB deacylase Mutant 1 5 with respect to parental sequences 
of Mutant 7-2 and Mutant 16. 

Figures 13 A, B, C, D and E show the DNA sequence of A. uiahensis ECB 
deacylase gene mutant M-l 5 genes created by in vitro heteroduplex formation followed 
by in vivo repair (SEQ ID NO:7). 
15 Figure 14 illustrates the process used for Example 3 to recombine 

mutations in RC1 and RC2 to yield thermostable subtilisin E. 

Figure 15 illustrates the sequences of RC1 and RC2 and the ten clones 
picked randomly from the transformants of the reaction products oi duplex formation as 
described in Example 3. The x's correspond to base positions that differ between RC1 and 
20 RC2. The mutation at 995 corresponds to amino acid substitution at 181, while that at 
1 107 corresponds to an ammo acid substitution at 218 in the subtilisin protein sequence. 

Figure 16 shows the results of screening 400 clones from the library 
created by heteroduplex formation and repair for initial activity (AO and residual activity 
(A r ). The ratio Aj/ A r was used to estimate the enzymes' thermostability. Data from 
25 active variants are sorted and plotted in descending order. Approximately 12.9% of the 
clones exhibit a phenotype corresponding to the double mutant containing both the 
N 1 8 1 D and the N2 1 8S mutations. 



DEFINITIONS 

30 Screening is. in general, a two-step process in which one first physically 

separates the cells and then determines which cells do and do not possess a desired 
property. Selection is a form of screening in which identification and physical separation 
are achieved simultaneously by expression of a selection marker, which, in some genetic 
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circumstances, allows ceils expressing the marker to survive while other cells die tor vice 
versa']. Exemplary screening members include iuciferase. 5galactosidase and green 
fluorescent protein. Selection markers include drug and toxin resistance genes. Although 
spontaneous selection can and does occur in the course of natural evolution, in the present 
methods selection is performed by man. 

An exogenous DNA segment is one foreign (or heterologous) to the cell or 
homologous to the cell but in a position within the host cell nucleic acid in which the 
element is not ordinarily found. Exogenous DNA segments are expressed to yield 
exogenous polypeptides. 

The term eene is used broadly to refer to any segment of DNA associated 
with a biological function. Thus, genes include coding sequences and/or the regulatory 
sequences required for their expression. Genes also include nonexpressed DNA segments 
that, for example, form recognition sequences for other proteins. 

The term "wild-type" means that the nucleic acid fragment does not 
> comprise any mutations. A "wild-type" protein means that the protein will be active at a 
level of activity found in nature and typically will comprise the amino acid sequence 
found in nature. In an aspect, the term "wild type" or "parental sequence" can indicate a 
starting or reference sequence prior to a manipulation of the invention. 

"Substantially pure" means an object species is the predominant species 
0 present (i.e.. on a molar basis it is more abundant than any other individual 

macromolecular species in the composition ), and preferably a substantially purified 
fraction is a composition wherein the object species comprises at least about 50 percent 
(on a molar basis) of all macromolecular species present. Generally, a substantially pure 
composition will compnse more than about 80 to 90 percent of all macromolecular 
.5 species present in the composition. Most preferably, the object species is purified to 
essential homogeneity (contaminant species cannot be detected in the composition by 
conventional detection methods ) wherein the composition consists essentially of a single 
macromolecular species. Solvent species, small molecules (<500 Daltons). and elemental 
ion species are not considered macromolecular species. 
30 Percentage sequence identity is calculated by comparing two optimally 

aliened sequences over the window of comparison, determining the numocr of positions 
at which the identical nucleic acid base occurs in both sequences to yield the number of 
matched positions, dividing the number of matched positions by the total number of 
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positions in the window of comparison. Optimal alignment of sequences for aligning a 
comparison window can be conducted by computerized implementations of algorithms 
GAP, BESTFIT, FASTA. and TFASTA in the Wisconsin Genetics Software Package 
Release 7.0, Genetics Computer Group, 575 Science Dr., Madison. WI. 

The term naturally-occurring is used to describe an object that can be 
found in nature as distinct from being artificially produced by man. For example, a 
polypeptide or polynucleotide sequence that is present in an organism (including viruses) 
that can be isolated from a source in nature and which has not been intentionally modified 
by man in the laboratory is naturally-occurring. Generally, the term naturally-occurring 
refers to an object as present in a non-pathological (undiseased) individual, such as would 
be typical for the species. 

A nucleic acid is operably linked when it is placed into a functional 
relationship with another nucleic acid sequence. For instance, a promoter or enhancer is 
operably linked to a coding sequence if it increases the transcription of the coding 
5 sequence. Operably linked means that the DNA sequences being linked are typically 
contiguous and, where necessary to join two protein coding regions, contiguous and in 
reading frame. However, since enhancers generally function when separated from the 
promoter by several kilobases and intronic sequences may be of variable lengths, some 
polynucleotide elements may be operably linked but not contiguous. 
0 A specific binding affinity between, for example, a ligand and a receptor, 

means a binding affinity of at least 1 x 10 6 M" 1 . 

The term "cognate" as used herein refers to a gene sequence that is 
evolutionarily and functionally related between species. For example but not limitation, 
in the human genome, the human CD4 gene is the cognate gene to the mouse CD4 gene, 
25 since the sequences and structures of these two genes indicate that they are highly 
homologous and both genes encode a protein which functions in signaling T cell 
activation through MHC class II-restricted antigen recognition. 

The term ■'heteroduplex* , refers to hybrid DNA generated by base pairing 
between complementary single strands derived from the different parental dupiex 
30 molecules, whereas the term "homoduplex" refers to double-stranded DNA generated b 
base pairing between complementary single strands derived from the same parental 
duplex molecules. 
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acquisition of a desired property. The substrates for the method are a populationof a, 
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II. Substrates For Shuffling 

The substrates for shuffling are variants of a reference polynucleotide that 
show some region(s) of similarity with the reference and other region(s) or point(s) of 
divergence. Regions of similarity should be sufficient to support annealing of 
polynucleotides such that stable heteroduplexes can be formed. Variants forms often 
show substantial sequence identity with each other (e.g., at least 50%. 75%. 90% or 99%). 
There should be at least sufficient diversity between substrates that recombination can 
generate more diverse products than there are starting materials. Thus, there must be at 
least two substrates differing in at least two positions. The degree of diversity depends on 
the length of the substrate being recombined and the extent of the functional change to be 
evolved. Diversity at between 0.1-25% of positions is typical. Recombination of 
mutations from very closely related genes or even whole sections of sequences from more 
distantly related genes or sets of genes can enhance the rate of evolution and the 
acquisition of desirable new properties. Recombination to create chimeric or mosaic 
genes can be useful in order to combine desirable features of two or more parents into a 
single gene or set of genes, or to create novel functional features not found in the parents. 
The number of different substrates to be combined can vary widely in size from two to 
10, 100, 1000, to more than 10 5 , 10 7, or 10 9 members. 

The initial small population of the specific nucleic acid sequences having 
mutations may be created by a number of different methods. Mutations may be created 
by error-prone PCR. Error-prone PCR uses low-fidelity polymerization conditions to 
introduce a low level of point mutations randomly over a long sequence. Alternatively, 
mutations can be introduced into the template polynucleotide by oligonucleotide-directed 
5 mutagenesis. In oligonucleotide-directed mutagenesis, a short sequence of the 

polynucleotide is removed from the polynucleotide using restriction enzyme digestion 
and is replaced with a synthetic polynucleotide in which various bases have been altered 
from the original sequence. The polynucleotide sequence can also be altered by chemical 
mutagenesis. Chemical mutagens include, for example, sodium bisulfite, nitrous acid. 
0 hvdroxvlarmne. hydrazine or formic acid. Other agents which are analogues of 
nucleotide precursors include nitrosoguanidine. 5-bromouraciL 2-aminopurine. or 
acridine. Generally, these agents are added to the PCR reaction in place of the nucleotide 
precursor thereby mutating the sequence. Intercalating agents such as proflavine. 
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acnfiavme. qumacnne and the like can also be used. Random mutagenesis of the 
polynucleotide sequence can also be achieved by irradiation with X-rays or ultraviolet 
lieht. Generally, plasmid DNA or DNA fragments so mutagenized are introduced into E. 
coli and propagated as a pool or library of mutant plasmids. 
5 Alternatively the small mixed population of specific nucleic acids can be 

found in nature in the form of different alleles of the same gene or the same gene from 
different related species (i.e.. cognate genes). Alternatively, substrates can be related but 
nonallelic genes, such as the immunoglobulin genes. Diversity can also be the result of 
previous recombination or shuffling. Diversity can also result from resynthesizmg genes 
1 0 encoding natural proteins with alternative codon usage. 

The starting substrates encode variant forms of sequences to be evolved. 
In some methods, the substrates encode variant forms of a protein for which evolution of 
a new or modified property is desired. In other methods, the substrates can encode 
variant forms of a plurality of genes constituting a multigene pathway. In such methods. 
15 variation can occur in one or any number of the component genes. In other methods, 
substrates can contain variants segments to be evolved as DNA or RNA binding 
sequences. In methods, in which starting substrates containing coding sequences, any 
essential regulatory sequences, such as a promoter and polyadenylation sequence, 
required for expression may also be present as a component of the substrate. 
20 Alternatively, such regulatory sequences can be provided as components of vectors used 
for cloning the substrates. 

The starting substrates can vary in length from about 50, 250. 1000, 
10,000, 100,000. 10 6 or more bases. The starting substrates can be provided in double- or 
sinele-stranded form. The starting substrates can be DNA or RNA and analogs thereof. 
If DNA. the starting substrates can be genomic or cDNA. If the substrates are RNA. the 
substrates are tvpically reverse-transcribed to cDNA before heteroduplex formation. 
Substrates can be provided as cloned fragments, chemically synthesized fragments or 
PCR amplification products. Substrates can derive from chromosomal, plasmid or viral 
sources. In some methods, substrates are provided in concatemeric form. 



J.5 



III. Procedures for Generating Heteroduplexes 

Heteroduplexes are generated from double stranded DNA substrates, by 
denaturing the DNA substrates and incubating under annealing conditions. Hybridization 
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conciitions for heteroduplex formation are sequence-dependent and are different in 
different circumstances. Longer sequences hybridize specifically at higher temperatures. 
Generally, hybridization conditions are selected to be about 25°C lower than the thermal 
melting point (Tm) for the specific sequence at a defined ionic strength and pH. The Tm 
5 is the temperature (under defined ionic strength, pH, and nucleic acid concentration) at 
which 50% of the probes complementary to the target sequence hybridize to the target 
sequence at equilibrium. 

Exemplary conditions for denaturation and renaturation of double stranded 
substrates are as follows. Equimolar concentrations (~ 1.0 - 5.0 nM) of the substrates are 
10 mixed in 1 x SSPE buffer (1 80 mM NaCL 1.0 mM EDTA, 10 mMNaH 2 P0 4 , pH 7.4) 
After heating at 96°C for 10 minutes, the reaction mixture is immediately cooled at 0°C 
for 5 minutes; The mixture is then incubated at 68°C for 2-6 hr. Denaturation and 
reannealing can also be carried out by the addition and removal of a denaturant such as 
NaOH. The process is the same for single stranded DNA substrates, except that the 
1 5 denaturing step may be omitted for short sequences. 

By appropriate design of substrates for heteroduplex formation, it is 
possible to achieve selection for heteroduplexes relative to reformed parental 
homoduplexes. Homoduplexes merely reconstruct parental substrates and effectively 
dilute recombinant products in subsequent screening steps. In general, selection is 
20 achieved by designing substrates such that heteroduplexes are formed in open-circles, 

whereas homoduplexes are formed as linear molecules. A subsequent transformation step 
results in substantial enrichment (e.g., 100-fold) for the circular heteroduplexes. 

Figure 1 shows a method in which two substrate sequences in separate 
vectors are PCR-amplified using two different sets of primers (PL P2 and P3. P4). 
25 Typically, first and second substrates are inserted into separate copies of the same vector. 
The two different pairs of primers initiate amplification at different points on the two 
vectors. Fig. 1 shows an arrangement in which the P1/P2 pnmer pairs initiates 
amplification at one of the two boundaries of the vector with the substrate and the P1/P2 
primer pair initiates replication at the other boundary in a second vector. The two primers 
30 in each primer pair prime amplification in opposite directions around a circular piasmid. 
The amplification products generated by this amplification are double-stranded linearized 
vector molecules in which the first and second substrates occur at opposite ends of the 
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vector. The amplification products are mixed, denatured and annealed. Mixing and 
denaturation can be performed in either order. Reanneaimg generates two linear 
homodupiexes. and an open circular heterodupiex containing one nick in each strand, at 
the initiation point of PCR amplification. Introduction of the amplification products into 
? host cells selects for the heteroduplexes relative to the homodupiexes because the former 
transform much more efficiently than the latter. 

It is not essential in the above scheme that amplification is initiated at the 
interface between substrate and the rest of the vector. Rather, amplification can be 
initiated at any points on two vectors bearing substrates provided that the amplification is 
0 initiated at different points between the vectors. In the general case, such amplification 
cenerates two linearized vectors in which the first and second substrates respectively 
occupy different positions relative to the remainder of the vector. Denaturation and 
reannealing generator heteroduplexes similar to that shown in Fig. 1. except that the nicks 
occur within the vector component rather than at the interface between plasmid and 
5 substrate. Initiation of amplification outside the substrate component of a vector has the 
advantage that it is not necessary to design primers specific for the substrate borne by the 
vector. 

Althouah Fig. 1 is exemplified for two substrates, the above scheme can be 
extended to any number of substrates. For example, an initial population of vector 
20 bearing substrates can be divided into two pools. One pooi is PCR-amplified from one 

set of primers, and the other pool from another. The amplification products are denatured 
and annealed as before. Heteroduplexes can form containing one strand from any 
substrate in the first pool and one strand from any substrate in the second pool. 
Alternatively, three or more substrates cloned into multiple copies of a vector can be 
25 subjected to amplification with amplification in each vector starting at a different point. 
For each substrate, this process generates amplification products varying in how flanking 
vector DNA is divided on the two sides of the substrate. For example, one amplification 
product might have most of the vector on one side of the substrate, another amplification 
product might have most of the vector on the other side of the substrate, and a further 
30 amplification product might have an equal division of vector sequence Hanking the 
substrate. In the subsequent annealing step, a strand of substrate can form a circular 
heterodupiex with a strand of any other substrate, but strands of the same substrate can 
oniv reanneal with each other to form a linear homodupiex. In a still further variation. 
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multipie substrates can be performed by performing multiple iterations of the scheme in 
Fig. 1 . After the first iteration, recombinant polynucleotides in a vector, undergo 
heteroduplex formation with a third substrate incorporated into a further copy of the 
vector. The vector bearing the recombinant polynucleotides and the vector bearing the 
third substrate are separately PCR amplified from different primer pairs. The 
amplification products are then denatured and annealed. The process can be repeated 
further times to allow recombination with further substrates. 

An alternative scheme for heteroduplex formation is shown in Fig. 2. 
Here, first and second substrates are incorporated into separate copies of a vector. The 
two copies are then respectively digested with different restriction enzymes. Fig. 2 shows 
an arrangement in which, the restriction enzymes cut at opposite boundaries between 
substrates and vector, but all that is necessary is to use two different restriction enzymes 
that cut at different places. Digestion generates linearized first and second vector bearing 
first and second substrates, the first and second substrates occupying different positions 
relative to the remaining vector sequences. Denaturation and reannealing generates open 
circular heteroduplexes and linear homoduplexes. The scheme can be extended to 
recombination between more than two substrates using analogous strategies to those 
described with respect to Fig. 1 . In one variation, two pools of substrates are formed, and 
each is separately cloned into vector. The two pools are then cute with different enzymes, 
and annealing proceeds as for two substrates. In another variation, three or more 
substrates can be cloned into three or more copies of vector, and the three or more result 
molecules cut with three or more enzymes, cutting at three or more sites. This generates 
three different linearized vector forms differing in the division of vector sequences 
flanking the substrate moiety in the vectors. Alternatively, any number of substrates can 
be recombined pairwise in an iterative fashion with products of one round of 
recombination annealing with a fresh substrate in each round. 

In a further variation, heteroduplexes can be formed from substrates 
molecules in vector-free form, and the heteroduplexes subsequently cloned into vectors. 
Such can be achieved by asymmetric amplification of first and second substrates as 
shown in Fig. 3. Asymmetric or single primer PCR amplifies only one strand of a duplex. 
By appropriate selection of primers, opposite strands can be amplified from two different 
substrates. On reannealing amplification products, heteroduplexes are formed from 
opposite strands of the two substrates. Because only one strand is amplified from each 
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be formed by d.rec.ona, c,on,n g . Two subs.ra.es for he.eroduplex formauon can be 
„b,ained bv PCR ampl.ficanon of chromosomal DNA and.o.ned .o opposne enos of a 
„ n ear vecor. D.reconaUlon.ng can be ach.eved by d.gesrma .he vecor ™.h nvo 
-,„ fferen, enzvmes. and d, g es„n g or adap.,n g f.r, and second subs.ra.es ,o be respecuveK 
• c mpauble ,„h cohes.ve ends of onlv »■ ,he ,vvo enzy mes used ,o cu, .he vecor. Tne 
firs, and second suos.ra.es can ,h„s be „ g a,ed a, opposne ends of a hnearrzed vecor 
fr aame„, Th.s scheme can oe ex.ended ,0 any number o, suos.ra.es by us,n 6 pr.nco.es 
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analogous to those described for Fie. 1 . For example, substrates can be divided into two 
pools before ligation to the vector. Alternatively, recombinant products formed by 
heterodupiex formation of first and second substrates, can subsequently undergo 
heteroduplex formation with a third substrate. 

5 

IV. Vectors and Transformation 

In general, substrates are incorporated into vectors either before or after 
the heteroduplex formation step. A variety of cloning vectors typically used in genetic 
engineering are suitable. 
1 o The vectors containing the DN A segments of interest can be transferred 

into the host cell by standard methods, depending on the type of cellular host. For 
example, calcium chloride transformation is commonly utilized for prokaryotic cells, 
whereas calcium phosphate treatment. Lipofection, or electroporation may be used for 
other cellular hosts. Other methods used to transform mammalian cells include the use of 
15 Polybrene, protoplast fusion, liposomes, electroporation, and microinjection, and 

biolisitics (see, generally, Sambrook et al., supra). Viral vectors can also be packaged in 
vitro and introduced by infection. The choice of vector depends on the host cells. In 
general, a suitable vector has an origin of replication recognized in the desired host cell, a 
selection maker capable of being expressed in the intended host cells and/or regulatory 
20 sequences to support expression of genes within substrates being shuffled. 

V. Types of Host Cells 

In general any type of cells supporting DNA repair and replication of 
heteroduplexes introduced into the cells can be used. Cells of particular interest are the 

25 standard cell types commonly used in genetic engineering, such as bacteria, particularly, 
E. coli (16. 17). Suitable E. coli strains include E. coli mutS, mutL. dam', and/or recA* . 
E.coli XL-10-Gold {[Tef A(mcrA) 1 83 AimcrCB-hsdSMR-mrr) 1 7 3 endAl supE44 ihi-1 
rccAl gyrA96 relAl lac Hte] [F* proAB lacFZAMIS TnlO (Tct r ) Amy Cam]), E. coli 
ES1301 mutS [Genotype: iacZ53. mutS201 : :Tn5, ihyA36. rha-5. mciBL deoC, INtrrnD- 

30 rmE)] (20. 24. 28-42). Preferred E. coli strains are E.coli SCSI 1 0 [Genotype: rpsi 
(Sir), thr. leu, enda. ihi-1, lacy, galk. gait, ara tona, tsx. dam. dcm. supE44, Aflac- 
proAB), [F. traD36, proA' B~ lacVZASil 5\ which have normal cellular mismatch repair 
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systems iP i. This strain type repairs mismatches and unmatches in the heteroduplex 
with little strand-specific preference, runner, because this strain is dam and dcm . 
piasmid isolated from the strain is unmethylated and therefore particularly amenable for 
further rounds of DNA duplex formation-mismatch repair (see below). Other suitable 
5 bacterial cells include gram-negative and gram-positive, such as Bacillus. Pseudomonas. 
and Salmonella. 

Eukaryotic organisms are also able to cam' out mismatch repair (43-48). 
Mismatch repair systems in both prokaryotes and eukaryotes are thought to play an 
important role in the maintenance of genetic fidelity during DNA replication. Some of 
1 0 the genes that play important roles in mismatch repair in prokaryotes. particularly mutS 

and mutL. have homologs in eukaryotes. in the outcome of genetic recombinations, and in 
genome stability. Wild-type or mutant 5. cercvisiae has been shown to carry out 
mismatch repair of heteroduplexes (49-56), as have COS-1 monkey cells (57). Preferred 
strains of yeast are Picchia and Saccharomyces. Mammalian cells have been shown to 
15 have the capacity to repair G-T to G-C base pairs by a short-patch mechanism (38. 58- 

63). Mammalian cells (e.g., mouse, hamster, primate, human), both cell lines and primary 
cultures can also be used. Such cells include stem cells, including embryonic stem cells, 
zygotes, fibroblasts, lymphocytes. Chinese hamster ovary (CHO). mouse fibroblasts 
(NIH3T3). kidney, liver, muscle, and skin cells. Other eucaryotic cells of interest include 
20 plant cells, such as maize, rice, wheat, cotton, soybean, sugarcane, tobacco, and 

arabidopsis: fish, algae, fungi (aspergillus. podospora. neurospora). insect (e.g.. baculo 
lepidoptera) (see. Winnacker. "From Genes to Clones." VCH Publishers. N.Y.. ( 1987), 
which is incorporated herein by reference ). 

In vivo repair occurs in a wide variety of prokaryotic and eukaryotic cells. 
25 Use of mammalian cells is advantage in certain application in which substrates encode 

polypeptides that are expressed only in mammalian cells or which are intended for use in 
mammalian cells. However, bacterial and yeast cells are advantageous for screening 
large libraries due to the higher transformation frequencies attainable in these strains. 

30 V. In Vitro DNA Repair Svstems 

As an alternative to introducing annealed products into host cells, annealed 
products can be exposed a DNA repair system m vitro. The DNA repair system can be 
obtained as extracts from repair-competent E. coii. yeast or any other cells ,64-67). 
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Repair-competent cells are lysed in appropriate buffer and supplemented with 
nucleotides. DNA is incubated in this cell extract and transformed into competent cells 
for replication. 

5 VI. Screening and Selection 

After introduction of annealed products into host cells, the host cells are 
typically cultured to allow repair and replication to occur and optionally, for genes 
encoded by polynucleotides to be expressed. The recombinant polynucleotides can be 
subject to further rounds of recombination using the heteroduplex procedures described 
10 above, or other shuffling methods described below. However, whether after one cycle of 
recombination or several, recombinant polynucleotides are subjected to screening or 
selection for a desired property. In some instances, screening or selection in performed in 
the same host cells that are used for DNA repair. In other instances, recombinant 
polynucleotides, their expression products or secondary metabolites produced by the 
15 expression products are isolated from such cells and screened in vitro. In other instances, 
recombinant polynucleotides are isolated from the host cells in which recombination 
occurs and are screened or selected in other host cells. For example, in some methods, it 
is advantageous to allow DNA repair to occur in a bacterial host strain, but to screen an 
expression product of recombinant polynucleotides in eucaryotic cells. The recombinant 
20 polynucleotides surviving screening or selection are sometimes useful products in 

themselves. In other instances, such recombinant polynucleotides are subjected to further 
recombination with each other or other substrates. Such recombination can be effected by 
the heteroduplex methods described above or any other shuffling methods. Further 
round(s) of recombination are followed by further rounds of screening or selection on an 
25 iterative basis. Optionally, the stringency of selection can be increased at each round. 

The nature of screening or selection depends on the desired property 
sought to be acquired. Desirable properties of enzymes include high catalytic activity, 
capacity to confer resistance to drugs, high stability, the ability to accept a wider (or 
narrower) range of substrates, or the ability to function in nonnatural environments such 
30 as organic solvents. Other desirable properties of proteins include capacity to bind a 
selected target, secretion capacity, capacity to generate an immune response to a given 
tareet. lack of immunogenicity and toxicity to pathogenic microorganisms. Desirable 
properties of DNA or RNA polynucleotides sequences include capacity to specifically 
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bind a civen protein target, and capacity to reguiate expression of operably linked coding 
sequences. Some of the above properties, such as drug resistance, can be selected by 
plating cells on the drug. Other properties, such as the influence of a regulatory sequence 
on expression, can be screened by detecting appearance of the expression product of a 
5 reporter gene linked to the regulatory sequence. Other properties, such as capacity of an 
expressed protein to be secreted, can be screened by FACS™. using a labelled antibody to 
the protein. Other properties, such as immunogenicity or lack thereof, can be screened by 
isolating protein from individual cells or pools of cells, and analyzing the protein in vitro 
or in a laboratory animal. 

10 

VU. Variations 

1. Demethvlation 

Most cell types methylate DNA in some manner, with the pattern of 
methylation differing between ceils types. Sites of methylation include 5-methylcytosine 
15 (m-O. N4-methylcytosine (m 4 C) and N 6 -methyladerune (m 6 A), 5-hydroxymethylcytosine 
(hm 5 C) and 5-hydroxymethyluracil (hnrUV In E. coli. methylation is effected by Dam 
and Dcm enzymes. The methylase specified by the dam gene methylates the N6-position 
of the adenine residue in the sequence G.VTC. and the methylase specified by the dcm 
gene methylates the C5 -position of the internal cytosine residue in the sequence 
20 CCWGG. DNA from plants and mammal is often subject to CG methylation meaning 
that CG or CNG sequences are methylated. Possible effects of methylated on cellular 
repair are discussed by references 18-20. 

In some methods. DNA substrates for heteroduplex formation are at least 
partially demethylated on one or both strands, preferably the latter. Demethvlation of 
25 substrate DNA promotes efficient and random repair of the heteroduplexes. In 

heteroduplexes formed with one strand dam-methylated and one strand unmethylated. 
repair is biased to the unmethylated strand, with the methylated strand serving as the 
template for correction. If neither strand is methylated, mismatch repair occurrs. but 
showes insignificant strand preference '.23. 24). 

Demethvlation can be performed in a variety of ways, in some methods, 
substrate DNA is demethyiated by PCR-ampiification. In some instances. DNA 
demethvlation is accomplished in one of the PCR steps in the heteroduplex formation 
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procedures described above. In other methods, an additional PCR step is performed to 
effect demethylation. In other methods, demethyiation is effected by passaging substrate 
DNA through methylation deficient host cells (e.g. an E. coli dam'dcm strain). In other 
methods, substrate DNA is demethylated in vitro using a demethylating enzymes. 
Demethylated DNA is used for heteroduplex formation using the same procedures 
described above. Heteroduplexes are subsequently introduced into DNA-repair-proficient 
but restriction-enzyme-defective cells to prevent degradation of the unmethylated 
heteroduplexes. 

2. Sealing Nicks 

Several of the methods for heteroduplex formation described above result 
in circular heteroduplexes beanng nicks in each strand. These nicks can be sealed before 
introducing heteroduplexes into host cells. Sealing can be effected by treatment with 
DNA ligase under standard ligating conditions. Ligation forms a phosphodiester bond to 
link two adjacent bases separated by a nick in one strand of double helix of DNA. 
Sealing of nicks increases the frequency of recombination after introduction of 
heteroduplexes into host cells. 

3. Error Prone PCR Attendant To Amplification 

Several of the formats described above include a PCR amplification step. 
Optionally, such a step can be performed under mutagenic conditions to induce additional 
diversity between substrates. 

VIII. Other Shuffling Methods 
5 The methods of heteroduplex formation described above can be used in 

conjunction with other shuffling methods. For example, one can perform one cycle of 
heteroduplex shuffling, screening or selection, followed by a cycle of shuffling by another 
method, followed by a further cycle of screening or selection. Other shuffling formats are 
described by WO 95/22625; US 5,605.793; US 5.81 1.238; WO 96/19256; Stemmer. 
0 Science 270. 1510 (1995); Stemmer et al.. Gene, 164. 49-53 (1995); Stemmer. 

Bio/Technology^ 13. 549-553 (1995); Stemmer. Proc. Natl. Acad. Sci. USA 91. 10747- 
10751 (1994); Stemmer. Nature 370. 389-391 (1994); Cramen et al.. Nature Medicine, 
2(1): 1-3. (1996); Cramen et aL. Nature Biotechnology 14. 315-319 ( 1 996); WO 98/42727; 
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WO 98/41622; WO 98/05764 and WO 98/42 7 2S. WO 98/27230 (each of which is 
incorporated by reference in its entirety for ail purposes). 

IX. Protein Analogs 

5 Proteins isolated by the methods also serve as lead compounds for the 

development of derivative compounds. The derivative compounds can include chemical 
modifications of amino acids or replace amino acids with chemical structures. The 
analogs should have a stabilized electronic configuration and molecular conformation that 
allows key functional groups to be presented in substantially the same way as a lead 

10 protein. In particular, the non-peptic compounds have spatial electronic properties which 
are comparable to the polypeptide binding region, but will typically be much smaller 
molecules than the polypeptides, frequently having a molecular weight below about 2 
CHD and preferably below about 1 CHD. Identification of such non-peptic compounds 
can be performed through several standard methods such as self-consistent field (CSF) 

1 5 analysis, configuration interaction (CHI) analysis, and normal mode dynamics analysis. 
Computer programs for implementing these techniques are readily available. See Rein et 
aU Computer-Assisted Modeling of Receptor -Ligand Interactions (Alan Liss, New York. 
1989). 



20 IX. Pharmaceutical Compositions 

Polynucleotides, their expression products, and secondary metabolites 
whose formation is catalyzed by expression products, generated by the above methods are 
optionally formulated as pharmaceutical compositions. Such a composition comprises 
one or more active agents, and a pharmaceutically acceptable earner. A variety of 
25 aqueous carriers can be used, e.g., water, buffered water, phosphate-buffered saline 

(PBS). 0.4% saline, 0.3% glycine, human albumin solution and the like. These solutions 
are sterile and generally free of particulate matter. The compositions may contain 
pharmaceutically acceptable auxiliary substances as required to approximate physio- 
locical conditions such as pH adjusting and buffering agents, toxicity adjusting agents 
30 and the like, for example, sodium acetate, sodium chloride, potassium chloride, calcium 
chloride and sodium is selected primarily based on fluid volumes, viscosities, and so 
forth, in accordance with the particular mode of administration selected. 
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EXAMPLES 

EXAMPLE 1 . Novel Rhizohium Flag Genes From Recombination Of Rhizobium Lupini 
Flag And Rhizobium Mcliloti FlaA 

Bacterial flagella have a helical filament, a proximal hook and a basal 
5 body with the flagellar motor (68). This basic design has been extensively examined in 
E. coli and S. typhimurium and is broadly applicable to many other bacteria as well as 
some archaea. The long helical filaments are polymers assembled from flageliin subunits, 
whose molecular weights range between 20,000 and 65,000. depending on the bacterial 
species (69). Two types of flagellar filaments, named plain and complex, have been 
10 distinguished by their electron microscopically determined surface structures (70). Plain 
filaments have a smooth surface with faint helical lines, whereas complex filaments 
exhibit a conspicuous helical pattern of alternating ridges and grooves. These 
characteristics of complex flagellar filaments are considered to be responsible for the 
brittle and (by implication) rigid structure that enables them to propel bacteria efficiently 
15 in viscous media (71-73). Whereas flagella with plain filaments can alternate between 

clockwise and counter clockwise rotation (68). all known flagella with complex filaments 
rotate only clockwise with intermittent stops (74). Since this latter navigation pattern is 
found throughout bacteria and archaea. it has been suggested that complex flagella may 
reflect the common background of an ancient, basic motility design (69). 
20 Differing from plain bacterial flagella in the fine structure of their 

filaments dominated by conspicuous helical bands and in their fragility, the filaments are 
also resistant against heat decomposition (72). Schmitt et al. (75) showed that 
bacteriophage 7-7-1 specifically adsorbs to the complex flagella of R. lupini H 13-3 and 
requires motility for a productive infection of its host. Though the flagellins from R. 
25 meliloti and R. lupini are quite similar, bacteriophage 7-7-1 does not infect R.meliloti. 
Until now complex flagella have been observed in only three species of soil bacteria: 
Pscudomonas rhodos (73), R.meliloti (76), and R.lupini HI 3-3 (70. 72). Cells of RAupini 
HI 3-3 posses 5 to 10 peritrichously inserted complex flagella. which were first isolated 
and analyzed by high resolution electron microscopy and by optical diffraction (70). 
30 Maruvama et al. (77) further found that a higher content of hydrophobic 

amino acid residues in the complex filament may be one of the main reasons for the 
unusual properties of complex flagella. By measuring mass per unit length and obtaining 
three-dimensional reconstruction from electron micrographs. Trachtenberc et al. (73. 78) 
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transformam library), (c) the sequence of ES01 (clone #1 from E. coli ESI 301 
iransformant library), and (d) the sequence of ES02 (clone #2 from E. coli ESI 301 
transformant library)- All four sequences were different from wild-type R. lupiniflaA and 
R. melilotiflaA sequences. Clones SCS02. ES01 and ES02 all contain a complete open- 
reading frame, but SCS01 was truncated. Figure 8 shows that recombination mainly 
occurred in the loop regions (unmatched regions). The flaA mutant library' generated 
from R. meliloti flaA and R. lupiniflaA can be transformed into E. coli SCSI 10, ESI 301, 
XLIO-Gold and JM109, and transformants screened for functional FlaA recombinants. 

EXAMPLE 2. Directed Evolution Of ECB Deacvlase For Variants With Enhanced 
Specific Activity 

Strepiomvccs are among the most important industrial microorganisms due 
to their ability to produce numerous important secondary metabolites (including many 
antibiotics) as well as large amounts of enzymes. The approach described here can be 
used with little modification for directed evolution of native Streptomyces enzymes, 
some or all of the genes in a metabolic pathways, as well as other heterologous enzymes 
expressed in Streptomyces. 

New antifungal agents are critically needed by the large and growing 
numbers of immune-compromised AIDS, organ transplant and cancer chemotherapy 
patients who suffer opportunistic infections. Echinocandin B (ECB), a lipopeptide 
produced by some species of Aspergillus, has been studied extensively as a potential 
antifungal. Various antifungal agents with significantly reduced toxicity have been 
generated by replacing the linoleic acid side chain of A. nidulans echinocandin B with 
different aryl side chains (79-83). The cyclic hexapeptide ECB nucleus precursor for the 
5 chemical acylation is obtained by enzymatic hydrolysis of ECB using Actinoplanes 

utahensis ECB deacvlase. To maximize the conversion of ECB into intact nucleus, this 
reaction is earned out at pH 5.5 with a small amount of miscible organic solvent to 
solubilize the ECB substrate. The product cyclic hexapeptide nucleus is unstable at pH 
above 5.5 during the long incubation required to fully deacylate ECB (84). The pH 
0 optimum of ECB deacvlase. however, is 8.0-8.5 and its activity is reduced at pH 5.5 and 
in the presence of more than 2.5% ethanol (84). To improve production of ECB nucleus 
it is necessary to increase the activity of the ECB deacvlase under these process-relevant 
conditions. 



WO 99 '29902 



PCT LS98 25698 



- 25 - 

Relatively iinle is known about ECB deacyiase. The enzyme is a 
heterociimer whose two subunits are cienved by processing of a singie precursor protein 
(83). The 19.9 kD a-subunh is separated from the 60.4 kD D-subumt by a 15-ammo acid 
spacer peptide that is removed along with a signal peptide and another spacer peptide in 
5 the native organism. The polypeptide is also expressed and processed into functional 

enzyme in Streptomyces lividans. the organism used for large-scale conversion of ECB by 
recombinant ECB deacyiase. The three-dimensional structure of the enzyme has not been 
determined, and its sequence shows so little similarity to other possibly related enzymes 
such as penicillin acylase that a structural model reliable enough to guide a rational effort 
0 to engineer the ECB deacyiase will be difficult to build. We therefore decided to use 
directed evolution (85) to improve this important activity. 

Protocols suitable for mutagenic PCR and random-priming recombination 
of the 2.4 kb ECB deacyiase gene (73% G-rC) have been descnbed recently (86). Here, 
we further describe the use of heteroduplex recombination to generate new ECB 
5 deacyiase with enhanced specific activity. 

In this case, two Actinoplanes utahensis ECB deacyiase mutants. M7-2 
and Ml 6, which show higher specific activity at pH 5.5 and in the presence of 10% 
MeOH were recombined using technique of the in vitro heteroduplex formation and in 
vivo mismatch repair . 

20 Figure 12 shows the physical maps of plasmids pM7-2 and pM16 which 

contain the genes for the M7-2 and Ml 6 ECB deacyiase mutants. Mutant M7-2 was 

obtained through mutagenic PCR performed directly on whole Streptomyces lividans cells 

containing wild-type ECB deacyiase gene, expressed from plasmid pSHP 150-2*. 

Streptomyces with pM7-2 show 1.5 times the specific activity of cells expressing the 
25 wild-type ECB deacyiase (86). Clone pM16 was obtained using the random-pnmmg 

recombination technique as descnbed (86. 87). It shows 2.4 times specific activity of the 

wild-type ECB deacyiase clone. 

A. Methods : 

, Q M7 _2 and M16 piasmid DNA (p.\r-2 and pM16x'Fie. 9) were purified 

from E. coii SCS110 (in separate reactions i. About 5.0 ug of unmethyiated M7-2 and 
M16 DNA were digested with A7jo 1 and Psh Al. respectively, at 37°C for 1 hour (Fig. 
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10). After agarose gel separation, me linearized DNA was purified using a Wizard PCR 
Prep Kit (Promega, WI. USA). 

Equimolar concentrations (2.0 nM) of the linearized unmethylated pM7-2 and pM16 
DNA were mixed in 1 x SSPE buffer (Ix SSPE: 180rnM NaCL 1.0 mM EDTA, lOmM 
5 NaH 2 P0 4 , pH 7.4). After heating at 96 °C for 10 minutes, the reaction mixture is 
immediately cooled at 0 °C for 5 minutes. The mixture was incubated at 68 °C for 3 
hours to promote formation of heteroduplexes. 

One microliter of the reaction mixture was used to transform 50 fil of 
E.coli ESI 301 mutS. SCSI 10 and JM109 competent cells. All transformants from E. coli 
10 ESI 301 mulS were pooled and E. coli SCSI 10 were pooled, A plasmid pool was isolated 
from each pooled library, and this pool was used to transform S. lividans TK23 
protoplasts to form a mutant library for deacylase activity screening. 
Transformants from the S. lividans TK23 libraries were screened for ECB deacylase 
activity with an in situ plate assay. Transformed protoplasts were allowed to regenerate 
15 on R2YE agar plates for 24 hr at 30°C and to develop in the presence of thiostrepton for 
48 hours. When the colonies grew to the proper size, 6 ml of 0.7% agarose solution 
containing 0.5 mg/ml ECB in 0.1 M sodium acetate buffer (pH 5.5) was poured on top of 
each R2YE-agar plate and allowed to develop for 1 8-24 hr at 30°C. Colonies surrounded 
by a clearing zone larger than that of a control colony containing wild-type plasmid 
20 pSHP150-2*, were selected for further characterization. 

Selected transformants were inoculated into 20 ml medium containing 
thiostrepton and grown aerobically at 30°C for 48 hours, at which point they were 
analyzed for ECB deacylase activity using HPLC. 100 \x\ of whole broth was used for a 
reaction at 30 °C for 30 minutes in 0.1 M NaAc buffer (pH 5.5) containing 10% (v/v) 
25 MeOH and 200 jag/ml of ECB substrate. The reactions were stopped by adding 2.5 

volumes of methanol, and 20 p.1 of each sample were analyzed by HPLC on a 100 x 4.6 
mm polyhydroxyethyl aspartamide column (PolyLC Inc.. Columbia. MD, USA) at room 
temperature using a linear acetonitriie gradient starting with 50:50 of A:B (A = 93% 
acetonitrile. 0.1% phosphoric acid; B = 70% acetonitriie. 0.1% phosphoric acid) and 
30 ending with 30:70 of A:B in 22 mm at a flow rate of 2.2 ml/mm. The areas of the ECB 
and ECB nucleus peaks were calculated and subtracted from the areas of the 
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butter was _ ed m g column 

containing 1 .0 M NaCl. A linear grad.ent trom 0 to 1 .0 M rs w 

f n 7 m , /min The ECB deacvlase fraction eluting at 0., M 

stain, and the concentration was determined using 

(Hercules. CA. USA>^ ^ ^ _ ^ ^ d _ ^ ^ of lhe ECB 
FCB substrate (84V Four US of each punfied ECB deacyiase 

5.5> comam.na .0% <«v, MeOH and d.fferen, concentrations of ECB subs, a t. 
Lavs were petfood ,„ duphcare. The reacuons were stopped by add.ng 2.5 volumes 
^ lanol. Id ,he HP L C assays were ca m ed out as descrtbed above. Tbe absorbance 
were recorded, and ,he ,n,„a, rates were emulated by leas, sauares re S ress,on of 
Ih e „me progress curves from which ,h= Km and the Kca, were calculated. 

Ac „v,„es as a funcuon of pH were measured for ,h= punfied ECB 

H , ^. c <; c an d 6 (0 1 M acetate buffer); /. 7.x 8 
deacvlases at 30°C at different pH values. >. 5.3 and (U. 

• «•„■»■ Q and 10 (0 1 M carbonate buffer) using the HPLL 

and8.5(0.lMphospnatebuffer,.9andl0(0.1 

- j cm deacvlases were were determined at 30 LinU.i m 
assav Stabilities ol punned ECB deacvlases wc 

:L,c buffer ,pH 5 .5 ■ conra.n.n, ,0* memanol. Samples were w„hd„w. a, d^ 
ume .mervals. and ,he res.dual acuv.ty was measured in the same buf.er wuh ,h= HPLC 

assav described above. 
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B. Results 

Fig. 1 1 shows that after one round of applying this heteroduplex repair 
technique on the mutant M7-2 and Ml 6 genes, one mutant (Ml 5) from about 500 original 
transformants was found to possess 3.1 times the specific activity of wild-type. 
5 Wild type and evolved Ml 5 ECB deacylases were purified and their kinetic parameters 
for deacylation of ECB were determined by HPLC. The evolved deacylases Ml 5 has an 
increased catalytic rate constant, kc at by 205%. The catalytic efficiency (kc al /K m ) of M20 
is enhanced by a factor of 2.9 over the wild-type enzyme. 

Initial rates of deacylation with the wild type and Ml 5 at different pH 
10 values from 5 to 10 were determined at 200 ^g/ml of ECB. The recombined Ml 5 is 
more active than wild type at pH 5-8. Although the pH dependence of the enzyme 
activity in this assay is not strong, there is a definite shift of 1 .0-1.5 units in the optimum 
to lower pH, as compared to wild type. 

The time courses of deactivation of the purified ECB deacylase mutant 
15 Ml 5 was measured in 0.1 M NaAc (pH 5.5) at 30°C. No significant difference in 
stability was observed between wild type and mutant Ml 5. 

The DNA mutations with respect to the wild type ECB deacylase sequence 
and the positions of the amino acid substitutions in the evolved variants M7-2, M16 and 
Ml 5 are summarized in Figure 12. 

20 The heteroduplex recombination technique can recombine parent 

sequences to create novel progeny. Recombination of the M7-2 and M16 genes yielded 
Ml 5, whose activity is higher than any of its parents (Fid. 13). Of the six base 
substitutions in Ml 5, five (at positions a50. a 171. 357, pi 29 and £340) were inherited 
from M7-2, and the other one (p30) came from Ml 6. 

25 This approach provides an alternative to existing methods of DNA 

recombination and is particularly useful in recombining large genes or entire operons. 
This method can be used to create recombinant proteins to improve their properties or to 
study structure-function relationship. 
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fv AMPI.F. 3. NovH Thermo*" ". *nr,lius Subfile Snhtilisin E Variants 

Th ls example demonstrates the use w vara heteroduplex formation 
followed bv ,„ vivo repair for combining sequence information from two different 
s ~ a uences in order to improve the thermostability of Bacillus subtilis subtilisin E. 

Genes RC1 and RC2 encode thermostable B. sublilis subtilisin E variants 
(88, The mutations at base positions 1 107 in RC1 and 995 in RC2 (Figure 14). giving 
rise to amino acid substitutions Asn218/Ser (N218S) and Asnl81/Asp (N181 ID), lead to 
im provements in subtilisin E thermostability; the remaining mutations, both synonymous 
and nonsvnonvmous. have no detectable effects on thermostability. At 65°C. the single 
, 0 variants NI81D and N21 8S have approximately 3-fold and 2-fold longer half-lives. 

respective* than wild subtilisin E. and variants containing both mutations have half-lives 
that are 8-fold lon.er (88). The different half-live, in a population of subtilisin E variants 
can therefore be used to estimate the efficiency by which sequence information is 
combined. In particular, recombination between these two mutations (in the absence of 
, s point mutations affecting thermostability) should generate a library in wh.ch 25% of the 
population exhibits the thermos/ability of the double mutant. Similarly. 25% of the 
population should exhibit wild-type like stability, as Nl 81D and N21 8S are eliminated at 
equal frequency. We used the fractions of the recombined population as a diagnostic 

-in A. Methods 

The stratecv undertone this example is shown in Fig. 1 5 . 
Subtilisin E thermostable mutant genes RC1 and RC2 (Fig. 14) are 986-bp 
foments including 45 nt of subtilisin E prosequence. the enure mature sequence and 1 1 3 
nt after the stop codon. The genes were cloned between Bam HI and Nae 1 in E coWB. 
25 sum, shuttle vector P BE3. resulting in P BE3-1 and P BE3-2. respectively. Plasmid 
DN A pBE3-l and pBE3-2 was isolated from E.coli SCSI 10. 

About 5.0 of ummethylated P BE3-1 andpBE3-2 DNA were digested 
with Bam HI and Nde 1. respectively, at 37»C for 1 hour. After agarose gel separation, 
eouimolar concentrations <2.0nM, of the linearized unmethylated pBE3-l and pBE,-2 
,0 were mixed in 1 x SSPE buffer < 180 mM NaCl. 1.0 mM EDTA. 10 mM NaH.PO, P H 

74, After heatine at 96=C for 10 m.nutes. the reaction mixture was immediately cooled 
, t 0 °C 'or < min The mixture was incubated at 6S'C for 2 hr for heteroduplexes to torn, 



WO 99/29902 



PCT/US98/25698 



-30- 

One microliter of the reaction mixture was used to transform 50 ul of E. 
coli ES 1301 mutS.E. coli SCS 1 1 0 and E. co//HB101 competent cells. 

The transformation efficiency with E. coli HB101 competent cells was 
about ten times higher than that of E. coli SCSI 10 and 15 times higher than that of E. coli 
5 ESI 301 mutS. But in all these cases, the transformation efficiencies were 10-250 times 
lower than that of the transformation with closed, covalent and circular control pUC19 
plasmids. 

Five clones from E. coli SCSI 10 mutant library and five from E. coli 
ESI 301 mutS library were randomly chosen, and plasmid DNA was isolated using a 
1 0 QI Aprep spin plasmid miniprep kit for further DNA sequencing analysis. 

About 2.000 random clones from E. coli HB10 1 mutant library were 
pooled and total plasmid DNA was isolated using a QIAGEN-100 column. 0.5-4.0 ug of 
the isolated plasmid was used to transform Bacillus subtil is DB428 as described 
previously (88). 

1 5 About 400 transformants from the Bacillus subtilis DB428 library were 

subjected to screening. Screening was performed using the assay described previously 
(88), on succinyl-Ala-Ala-Pro-Phe-p-nitroanilide. B. subtilis DB428 containing the 
plasmid library were grown on LB plates containing kanamycin (20 |ag/ml) plates. After 
18 hours at 37°C single colonies were picked into 96-well plates containing 200 \i\ 
20 SG/kanamycin medium per well. These plates were incubated with shaking at 37°C for 
24 hours to let the cells to grow to saturation. The cells were spun down, and the 
supernatants were sampled for the thermostability assay. 

Two replicates of 96-well assay plates were prepared for each growth plate 
by transferring 1 0 \i\ of supernatant into the replica plates. The subtilisin activities were 
25 then measured bv adding 100 ul of activity assay solution (0.2 mM succinyl-Ala-Ala-Pro- 
Phe-p-nitroanilide. 100 mM Tris-HCl. 10 mM CaCh, pH 8.0. 37°C). Reaction velocities 
were measured at 405 nm to over 1.0 mm in a ThermoMax microplate reader (Molecular 
Devices. Sunnyvale CA). Activity measured at room temperature was used to calculate 
the fraction of active clones (clones with activity less than 10% of that of wild type were 
30 scored as inactive). Initial activity (A,) was measured after incubating one assay plate at 
65°C for 10 minutes by immediately adding 100 ul of prewarmed (37°C) assay solution 
(0.2 mM succmyl-Aia-Aia-Pro-Phe-p-nitroaniiidc. 100 mM Tris-HCl. pH 8.0. 10 mM 
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CaCK pH 8.0) into each we:!. Residual activity ( Ar) was measured after 40 minute 
incubation. 

B. Results 

In vitro heterodupiex formation and in vivo repair was earned out as 
5 described above. Five clones from E. coli SCSI 10 mutnnt library and five from E. coil 
ESI 301 mutS libraries were selected at random and sequenced. Fig. 14 shows that four 
out of the ten clones were different from the parent genes. The frequency of occurrence 
of a particular point mutation from parent RC1 or RC2 in the resulting genes ranged from 
0% to 50%. and the ten point mutations in the heterodupiex have been repaired without 
10 strong strand-specific preference. 

Since none of the ten mutations locates within the dem sue. the mismatch 
repair appears generally done via the E. coli lone-patch mismatch repair systems. The 
system repairs different mismatches in a strand-specific mariner using the state of N6- 
methylation of adenine in GATC sequences as the major mechanism for determining the 
1 5 strand to be repaired. With heterodupiexes methylated at GATC sequences on only one 
DN A strand, repair was shown to be highly biased to the unmethylated strand, with the 
methylated strand serving as the template for correction. If neither strand was 
methylated, mismatch repair occurred, but showed little strand preference (23. 24). These 
results shows that it is preferable to demethylate the DNA to be recombined to promote 
20 efficient and random repair of the heterodupiexes. 

The rates of subtilisin E thermo-mactivation at 65°C were estimated by 
analyzing the 400 random clones from the Bacillus subtilis DB42S library. The 
thermostabilities obtained from one 96-well plate are shown in Figure 16. plotted in 
descendine order. About 12.9% of the clones exhibited thermostability comparable to the 
25 mutant with the NISID and N218S double mutations. Since this rate is only half of that 
expected for random recombination of these two markers, it indicates that the two 
mismatches at positions 995 and 1 107 within the heterodupiexes have been repaired with 
lower position randomness. 

Sequence analysis of the clone exhibiting the highest thermostability 
30 amonc the screened 400 transTormams Irom the E. coli SCSI 10 heterodupiex library- 
confirmed the presence of both NISID and N218S mutations. Among the 400 
transformants from the B.subiilis DB428 library that were screened, approximately 91% 
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of the clones expressed N181D- and/or N218S-type enz\'me stabilities, while about 8.0% 
of the transformants showed only wild-type subtilisin E stability. 

Less than 1.0% inactive clone was found, indicating that few new point 
mutations were introduced in the recombination process. This is consistent with the fact 
that no new point mutations were identified in the ten sequenced genes (Figure 14). 
While point mutations may provide useful diversity for some in vitro evolution 
applications, they can also be problematic for recombination of beneficial mutations, 
especially when the mutation rate is high. 

EXAMPLE 4. Optimizing Conditions For The Heteroduplex Recombination. 

We have found that the efficiency of heteroduplex recombination can 
differ considerably from gene to gene [17,57]. In this example, we investigate and 
optimize a variety of parameters that improve recombination efficiency. 
DNA substrates used in this example were site-directed mutants of green fluorescent 
5 protein from Aequorea victoria. The GFP mutants had a stop codon(s) introduced at 

different locations along the sequence that abolished their fluorescence. Fluorescent wild 
type protein could be only restored by recombination between two or more mutations. 
Fraction of fluorescent colonies was used as a measure of recombination efficiency. 

0 A. Methods 

About 2-4 jag of each parent piasmid was used in one recombination 
experiment. One parent piasmid was digested with Pst I endonuclease another parent 
with EcoKl. Linearized plasmids were mixed together and 20 x SSPE buffer was added to 
the final concentration lx (180 mM NaCL 1 mM EDTA, 10 mM NaH 2 P0 4 , pH 7.4). The 
15 reaction mixture was heated at 96°C for 4 minutes, immediately transferred on ice for 4 
minutes and the incubation was continued for 2 hours at 68°C. 

Target genes were amplified in a PCR reaction with primers corresponding 
to the vector sequence of pGFP piasmid. Forward primer: 5'- 
CCGACTGGAAAGCGGGCAGTG-3\ reverse primer 5*- 
30 CGGGGCTGGCTTAACTATGCGG-3". PCR products were mixed together and purified 
using Qiagen PCR purification kit. Purified products were mixed with 20 x SSPE buffer 
and hybridized as described above. Annealed products were precipitated with ethanol or 
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punned on Qiagen columns and digested with EcoRl and Psil enzymes. Digested 
products were iigated into Pstl and EcoRl digested pGFP vector. 

dUTP was added into PCR reaction at final concentrations 200 uM. 40uM. 
8 uM. 1 .6 uM, 0.32 faM. PCR reaction and subsequent cloning procedures were 

performed as described above. 

Recombinant piasmids were transformed into XL 10 E. coli strain by a 
modified chemical transformation method. Cells were plated on ampicillin containing LB 
agar plates and grown overnight at 37°C, followed by incubation at room temperature or 
at 4°C until fluorescence developed. 

B. Results . 

1 . Effect of ligation on recombination efficiency . 

Two experiments have been performed to test the effect of breaks in the 
DNA heteroduplex on the efficiency of recombination. In one expenment heteroduplex 
5 plasmid was treated with DNA ligase to close all existing single-strand breaks and was 
transformed in identical conditions as an unligated sample (see Table 1 ). The Iigated 
samples show up to 7-fold improvement in recombination efficiency over unligated 
samples. 

In another experiment. dUTP was added into PCR reaction to introduce 
0 additional breaks into DNA upon repair by uracyl N-glycosylase in the host cells. Table 
2 shows that dUMP incorporation significantly suppressed recombination, the extent of 
suppression increasing with increased dUTP concentration. 

2. Effect of plasmid size on the efficiency of 
>5 heteroduplex formation . 

Plasmid size was a significant factor affecting recombination efficiency. 
Two plasmids pGFP (3.3 kb) and a Bacillus shuttle vector pCTl (about 9 kb) were used 
in preparing circular heteroduplex-like piasmids following traditional heteroduplex 
protocol. For the purpose of this expenment (to study the effect of plasmid size on 
30 duplex formation), both parents had the same sequences. While pGFP formed about 30 
40% of circular plasmid. the shuttle vector yielded less than 1 0° o oi this form. 
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Increase in plasmid size decreases concentration of the ends in the vicinity 
of each and makes annealing of very long (>0.8 kb) ends that are single-stranded more 
difficult. This difficulty is avoided by the procedure shown in Fig. 3, in which 
heteroduplex formation occurs between substrates in vector-free form, and, 
5 heteroduplexes are subsequently inserted into a vector. 

3. Efficiency of Recombination vs. Distance Between Mutations 

A series of GFP variants was recombined pairwise to study the effect of 
distance between mutations on the efficiency of recombination. Parental genes were 
10 amplified by PCR. annealed and ligated back into pGFP vector. Heteroduplexes were 
transformed into XL10 E.coli strain. 

The first three columns in Table 3 show the results of three independent 
experiments and demonstrate the dependence of recombination efficiency on the distance 
between mutations. As expected recombination becomes less and less efficient for very 
1 5 close mutations. 

However, it is still remarkable that long-patch repair has been able to 
recombine mutations separated by only 27 bp. 

The last line in Table 3 represents recombination between one single and 
one double mutants. Wild type GFP could only be restored in the event of double 
20 crossover with each individual crossover occurring in the distance of 99 bp only, 

demonstrating the ability of this method to recombine multiple, closeiy-spaced mutations. 

4. Elimination Of The Parental Double Strands 

From Heteroduplex Preparations. 
25 Annealing of substrates in vector-free form offers size-advantages relative 

to annealing of substrates as components of vectors, but does not allow selection for 
heteroduplexes relative to homoduplexes simply by transformation into host. 
Asymmetric PCR reactions with only one primer for each parent seeded with appropriate 
amount of previously amplified and purified gene fragment were run for 100 cycles. 
30 ensuring a 100-fold excess of one strand over another. Products of these asymmetrical 
reactions were mixed and annealed together producing only a minor amount of 
nonrecombinant duplexes. The last column in Table 3 shows the recombination 
efficiency obtained from these enriched heteroduplexes. Comparison of the first three 
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columns with the fourth one demonstrates the improvement achieved by asymmetric 
synthesis of the paremai strands. 

While the foregoing invention has been described in some detail for 
purposes of clarity and understanding, it will be clear to one skilled in the an from a 
reading of this disclosure that various changes in form and detail can be made without 
departing from the true scope of the invention. All publications and patent documents 
cited in this application are incorporated by reference in their entirety for all purposes to 
the same extent as if each individual publication or patent document were so individually 
denoted. 
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WHAT IS CLAIMED IS: 

1 1 A method for evolving a polynucleotide toward acquisition of a 

2 desired property, comprising 

(a) incubating a population of parental polynucleotide variants under 

4 conditions to generate annealed polynucleotides comprising heteroduplexes; 

5 (b) exposing the heteroduplexes to a cellular DNA repair system to 

6 convert the heteroduplexes to parental polynucleotide variants or recombined 

7 polynucleotide variants; 

8 (c) screening or selecting the recombined polynucleotide variants for the 

9 desired property. 

1 2. The method of claim 1. wherein the heteroduplexes are exposed to 

2 the cellular DNA repair system in vitro. 

1 3. The method of claim 2, wherein the cellular DNA repair system 

2 comprises cellular extracts. 

1 4. The method of claim 1 , further comprising introducing the 

2 heteroduplexes into cells, whereby the heteroduplexes are exposed to the DNA repair 

3 system of the cells in vivo. 

1 5. The method of claim 4. wherein the annealed polynucleotides 

2 further comprise homoduplexes and the introducing step selects for transformed cells 

3 comprising the heteroduplexes relative to transformed cells comprising homoduplexes. 

1 6. The method of claim 4, wherein a first polynucleotide variant is 

2 provided as a component of a first vector, and a second polynucleotide variant is provided 

3 as a component of a second vector, and the method further comprises convening the first 

4 and second vectors to linearized forms in which the first and second polynucleotide 

5 variants occur at opposite ends, whereby m the incubating step single-stranded forms of 

6 the first imeanzed vector reanneai with each other to form linear first vector, single- 

- stranded forms of the second linearized vector reanneai with each other to form linear 
8 second vector, and single-stranded Imeanzed forms of the first and second vectors anneal 
q with each to form a circular heteroduplex bearing a nick in each strand, and the 
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introducing step selects for transformed cells comprising the circular 
heteroduplexes relative to the linear first and second vector. 

7. The method of claim 6. wherein the first and second vectors axe 
convened to linearized forms by PCR. 



1 8. The method of claim 6, wherein the first and second vectors are 

2 converted to linearized forms by digestion with first and second restriction enzymes. 

1 9. The method of claim k wherein the population of polynucleotide 

2 variants are provided in double stranded form, and the method further comprising 

3 converting the double stranded polynucleotides to single stranded polynucleotides before 

4 the annealing step. 

1 10. The method of claim 1 , wherein the converting step comprises: 

2 conducting asymmetric amplification of the first and second double 

3 stranded polynucleotide variants to amplify a first strand of the first polynucleotide 

4 variant, and a second strand of the second polynucleotide variant, whereby the first and 

5 second strands anneal in the incubating step to form a heteroduplex. 

1 11. The method of claim 1 0. wherein the first and second double- 

2 stranded polyncueltoides variants are provided in vector-free form, and the method 

3 further comprises incorporating the heteroduplex into a vector. 

1 12. The method of claim 4 wherein the population of polynucleotides 

2 comprises first and second polynucleotides provided in double stranded form, and the 

3 method further comprises incorporating the first and second polynucleotides as 

4 components of first and second vectors, whereby the first and second polynucleotides 

5 occupy opposite ends of the first and second vectors, whereby in the incubating step 

6 single-stranded forms of the first linearized vector reanneal with each other to form linear 

7 first vector, single-stranded forms of the second linearized vector reanneal with each 

8 other to form linear second vector, and single-stranded linearized forms of the first and 

9 second vectors anneal with each to form a circular heteroduplex bearing a nick in each 
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1 strand, and the introducing step selects for transformed cells comprises the 

2 circular heteroduplexes relative to the linear first and second vector. 

1 13. The method of claim 4. further comprising sealing nicks in the 

2 heteroduplexes to form covalently-ciosed circular heteroduplexes before the introducing 

3 step. 

\ 1 4 . The method of claim 1 1 . wherein the first and second 

2 polynucleotides are obained from chromosomal DNA.. 

1 15. The method of claim 1. further comprising repeating steps (a)-(c) 

2 whereby the incubating step in a subsequent cycle is performed on recombinant vanants 

3 from a previous cycle. 

! 1 6. The method of claim 1 , wherein the polynucleotide vanants encode 

2 a polypeptide. 

1 1 7. The method of claim 1 . wherein the population of polynucleotide 

2 variants comprises at least 20 variants. 

! 18. The method of claim 1 . wherein the population of polynucleotide 

2 variants are at least 1 0 kb in length. 

! 19. The method of claim 1. wherein the population of polynucleotide 

2 variants comprises natural variants. 

! 20. The method of claim 1 . wherein the population of polynucleotides 

2 comprises variants generated by mutagenic PCR. 

\ 21. The method of claim 1. wherein the population of polynucleotide 

2 vanants comrises variants generated by site directed mutagenesis. 

22. The method of claim 1 . wherein the cells are bacterial cells. 

23 The method of claim 1 . further comprising at least partially 
2 demethylating the population of variant polynucleotides. 
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1 24. The method of claim 23. whether the at least partially 

2 demethylating step is performed by PCR amplification of the population of variant 

3 polynucleotides. 

1 25. The method of claim 23. wherein the at least partially 

2 demethvlating step is performed by amplification of the population of variant 

3 polynucelotides in host cells. 

1 26. The method of claim 25, wherein the host cells are defective in a 

2 gene encoding a methylase enzyme. 

1 27. The method of claim 1, wherein the population of variant 

2 polynucleotide variants comprises at least 5 polynucleotides having at least 90% sequence 

3 identity with one another. 

1 28. The method of claim 1 , further comprising isolating a screened 

2 recombinant variant. 

1 29. The method of claim 28. further comprising expressing a screened 

2 recombinant variant to produce a recombinant protein. 

1 30. The method of claim 29 further comprising formulating the 

2 recombinant protein with a carrier to form a pharmaceutical composition. 

1 31. The method of claim 1, wherein the polynucleotide variants encode 

2 enzymes selected from the group consisting of proteases, lipases, amylases, cutinases, 

3 cellulases, amylases, oxidases, peroxidases and phytases. 

1 32. The method of claim 1. wherein the polynucleotide variants encode 

2 a polypeptide selected from the group consisting of insulin. ACTH. glucagon. 

3 somatostatin, somatotropin, thymosin, parathyroid hormone, pigmentary hormones. 

4 somatomedin, erthropoietin. luteinizing hormone, chorionic gonadotropin, hyperthalmic 

5 reieasinn factors, antidiuretic hormones, thyroid stimulating hormone, reiaxin. interferon. 

6 thrombopoietin (TPCM. and prolactin. 
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33. The method of claim 1 . wherein the polynucleotide van ants 
encode a plurality of enzymes forming a metabolic pathway. 

34. The method of claim L wherein the polynucleotide variants are in 
concatemeric form. 
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